Analytically solvable model of an electronic Mach-Zehnder interferometer 
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We consider a class of models of non-equilibrium electronic Mach-Zehnder interferometers built on 
integer quantum Hall edges states. The models are characterized by the electron-electron interaction 
being restricted to the inner part of the interferometer and transmission coefficients of the quantum 
quantum point contacts, defining the interferometer, which may take arbitrary values from zero 
to one. We establish an exact solution of these models in terms of single-particle quantities — 
determinants and resolvents of Fredholm integral operators. In the general situation, the results 
can be obtained numerically. In the case of strong charging interaction, the operators acquire the 
block Toeplitz form. Analyzing the corresponding Riemann-Hilbert problem, we reduce the result 
to certain singular single-channel determinants (which are a generalization of Toeplitz determinants 
with Fisher-Hartwig singularities), and obtain an analytic result for the interference current (and, 
in particular, for the visibility of Aharonov-Bohm oscillations). Our results, which are in good 
agreement with experimental observations, show an intimate connection between the observed "lobe" 
structure in the visibility of Aharonov-Bohm oscillations and multiple branches in the asymptotics 
of singular integral determinants. 

PACS numbers: 71.10.Pm, 73.23.-b, 73.43.-f, 85.35.Ds 



I. INTRODUCTION 

Electronic Mach-Zehnder interferometers (MZIs) real- 
ized with edge states in the integer quantum Hall (QH) 
regime have attracted a lot of attention recently because 
of a striking interplay between the quantum coherence 
and effects of electron-electron interaction observed in 
these mesoscopic devices^HE] By analogy to the op- 
tical interferometer, the chiral edge states in the elec- 
tronic MZI, playing the role of light beams, are coupled 
by quantum point contacts (QPCs), which act as elec- 
tron beam-splitters (see Fig. [l}. The differential conduc- 
tance measured in the above experiments shows strong 
Aharonov-Bohm (AB) oscillations. They are a manifes- 
tation of quantum coherence of electrons, propagating 
through different arms of interferometer, and are quan- 
tified in terms of visibility. The most remarkable exper- 
imental observation is that the out-of-equilibrium visi- 
bility does not decrease monotonically with voltage but 
rather demonstrates a sequence of decaying oscillations 
( "lobes" ) . Such a dependence cannot be explained within 
an assumption of non-interacting electrons. 

Investigation of quantum interference and decoherence 
in Aharonov-Bohm rings and interferometers has a long 
historj^. In particular, much attention has been paid 
to sources of dephasing that may arise from the exter- 
nal noise^H2H or are ^ ne resu it of the intrinsic electron- 
electron interactio n 22 " 26 !. The advent of QH interferome- 
ters has renewed the interest in this problem, with a con- 
siderable number of recent theoretical workg^EU aiming 
at a resolution of the "visibility puzzle" in MZIs. These 
recent theories can be subdivided int o the app roaches as- 
suming contact!^ 32 1 36 1 and long-rang d 27 ' 28133 " 3 ^ Coulomb 



interaction. Despite the fact that the model of contact 
e-e interaction may successfully describe the related ex- 
periments on the energy relaxation in the QH edge states 
at filling factor i/=sPtal results of Refs. [351 and H2 in- 
dicate that the account of the long-range character of 
Coulomb interaction is of central importance for a full 
understanding of non-equilibrium phenomena in MZIs. 

The natural choice of a theoretical approach to one- 
dimensional (ID) interacting electrons in the QH edge 
states is that of bosonizatiorP^J. However, in the case 
under interest, one faces serious complications when 
trying to apply this approach. First, already in the 
single-channel problem, the bosonized action of the 
theory ceases to be Gaussian under non-equilibrium 
conditions" (see also a related earlier worl^ on the 
non-equilibrium Fermi-edge singularity). This difficulty 
has been solved by development of the non-equilibrium 
bosonization formalism yielding results for physical ob- 
servables in terms of single-particle Fredholm determi- 
nants which are o f Toeplitz type for a short-range inter- 
action model 4445 . Second, an even more severe obstacle 
arises when one describes electron scattering at QPCs. 
Specifically, electron tunneling between two edge chan- 
nels yields the cos-like term in the bosonized Hamiltonian 
(an the action), impeding a solution to the problem. For 
this reason, almost all recent theories of MZIs consider 
the limit of weakly coupled edge states where the pertur- 
bative treatment of electron tunneling at QPCs is justi- 
fied. This is rather unfortunate since in the experiment 
transmission coefficients of both QPCs are usually close 
to one half. While in a very restricted set of models exact 
solutions via the Bethe ansatz are available^^, the sys- 
tems we are interested in do not belong to this class. On 



the other hand, it would be clearly highly advantageous 
to have an analytically treatable mode for integer QH 
MZI for an arbitrary number of edge channels, arbitrary 
interaction range and strength, and arbitrary transmis- 
sions at QPC. 

In this paper we consider the model of the MZI oper- 
ating at integer filling factor v where electrons interact 
only when they are inside the interferometer. The model 
is specified by two single-particle scattering matrices of 
the QPCs defining the interferometer and by the model of 
Coulomb interaction inside the interferometer. We focus 
on the model of "maximally long-range" interaction when 
the interaction energy depends only on total charges col- 
lected within each of the arms and is characterized by 
an electrostatic charging energy E c . Let us emphasize 
that this restriction is not crucial: within this approach 
one can, in principle, consider any interaction within the 
interior region of MZI. 

For the case v = 1, such a model was introduced in 
Refs. and 1301 In Ref- EU an exact solution to it at 
v = 1 was obtained by using a combination of bosoniza- 
tion and refermionization techniques and was expressed 
in terms of single-particle determinant and resolvent that 
were evaluated numerically. Our way to treat the prob- 
lem is different in many aspects. We consider a MZI 
with an arbitrary number of edge states and use the non- 
equilibrium functional bosonization approach developed 
by us previously^. Within this framework we demon- 
strate that an interfering current can be expressed in 
terms of a Fredholm functional determinant of the single- 
particle "counting" operator which bears resemblance to 
the problem of electron full counting statistics (FCS) 
of mesoscopic transportPl In general, this determinant 
should be evaluated numerically. In the limit of strong 
interaction E c 3> 1/r, where r is the electron flight 
time through the MZI, the "counting" operator takes the 
block Toeplitz form. Under this condition, a fully ana- 
lytical treatment turns out to be possible. By solving 
the Riemann-Hilbert problem, we get rid of the matrix 
structure and express the result in terms of a determinant 
of a single- channel singular integral operator generalizing 
Toeplitz determinants with Fisher-Hartwig singularities. 
Determinants of this type have been studied in Ref. 1521 
where a conjecture on their asymptotic behavior was for- 
mulated and supported by a large body of analytical and 
numerical arguments. This allows us to obtain the result 
in a closed analytic form. 

Our analytical result demonstrates that the "lobe" 
pattern in visibility is a many-body interference effect 
resulting from the quantum superposition of (at least) 
two many-particle scattering amplitudes with the mutual 
phase difference which is linear in voltage. In the limit 
of strong interaction we find the scaling exponents which 
describe power-law dependences of these amplitudes on 
voltage and obtain the non-equilibrium dephasing rate 
governing the exponential suppression of visibility with 
bias (or, equivalently, with the length of the arms of the 
MZI) . The power-law exponents as well as the dephasing 
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FIG. 1. Layout of an electronic Mach-Zehnder interferometer 
built on quantum-Hall edge states at filling factor v = 2. 
Quantum point contact (QPC1 and QPC2) characterized by 
transparencies Ti(2) are used to partially mix the outer edge 
channels. All Ohmic contacts are grounded, except for the 
source terminal S3 which is kept at voltage V. The current is 
measured in the drain terminal D\. The QPCO can be used to 
dilute the incoming current in the outer channel by changing 
the transparency To- 



rate depend on the transmission coefficient of the first 
QPC and on the filling factor v. 

Our analytical findings are corroborated and comple- 
mented by numerical evaluations of the Fredholm de- 
terminants determining the exact solution for arbitrary 
E c ~ 1/r. At E c 3> 1/t our numerical results pro- 
vide further support to the aforementioned conjecture 
of Ref. 52 At moderate charging energy E c > 1/r and 
v=1 the obtained results match very well experimental 
observations. 

The remainder of the paper is organized as follows. 
Sectionlnlis devoted to the exposition of our main results 



and their physical interpretation. In Sec. Ill we present 
the non-equilibrium functional bosonization approach. 
We demonstrate that the MZI problem defined above 
is exactly solvable by means of the instanton method. 



In the limit of strong interaction (Sec. IV I, the full an- 
alytical treatment becomes possible. It is based upon 
the asymptotic results for the generalized Toeplitz de- 
terminants. We show the relation of the MZI problem to 
the latter theory and evaluate analytically the Aharonov- 
Bohm conductance. In Sec. [V] we consider the influence 
of an additional quantum point contact (QPCO) diluting 
the incoming current in one of the channels and develop 
numerically exact approach to solve the problem in the 
case of a moderate charging energy. Finally, in Sec. |VI| wc 
summarize the findings of this work and discuss prospects 
for future research. 



II. RESULTS AND QUALITATIVE DISCUSSION 

In this section we set the stage by defining the theoret- 
ical model of the Mach-Zehnder interferometer and then 
present our results and give their physical interpretation. 
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This part of the paper is self-contained and can be read 
independently of other sections, where we provide tech- 
nical details of the calculations. 



A. Model of the MZI 

We consider the Mach-Zchndcr interferometer (MZI) 
realized with edge states in the QH regime at integer 
filling factor v. The experimental layout (in case of v = 
2) and the scheme of the MZI are shown in Figs. [I] and [2] 
In this setup the outer chiral channels, propagating along 
different arms of the MZI (we denote them by the index 
±) are coupled by means of two quantum point contacts 
(QPCs), located at points x ± . An additional quantum 
point contact, QPCO, is used to bias the incoming outer 
channel at the upper edge by voltage V (and in general 
also to dilute the incoming current if the transmission 
coefficient of the QPCO is tuned to a value T > 0). 
We further assume that all inner chiral channels are fully 
reflected from each QPC and, in particular, the incoming 
inner channels are grounded. This layout of the MZI and 
the bias scheme are realized in most of the experiments 
(an exception is Ref. [T^l where the MZI setup at v = 2 
did not contain additional QPCO). 

A theoretical model considered throughout the paper 
is specified by the action S = Sq + S- mt of interacting ID 
fermions, 



S = ^ dtdxip e (x)(id t + ivd x )ip e (x), (1) 
g=± J 



are inside the interferometer. Thus in Eq. ^ 

AA e = 2/ ^ kg (x + 0)^ kg (x) (3) 
k=i Jx l 

is the total number of electrons in the upper/lower arm of 
the MZI. The action A- ln t describes intra- and interchan- 
nel Coulomb interaction which is maximally non-local (or 
long-range) in space. At the same time the inter-edge 
interaction is disregarded in our model, which is moti- 
vated by the fact that different edges are spatially well 
separated. At the QPCs outgoing fermion fields (with 
channel index k = 1) are related to the incoming ones by 
the scattering matrices 



0) 



,1/2 rpl/2 \ 



-,1/2 



T 



iR 



1/2 



(4) 
(5) 



where Tj and Rj 



are reflection and transmission coeffi- 
cients at the j-th QPC. 

The model of the MZI with the above action A is ex- 
actly solvable for any value of the charging energy E c and 



transmission coefficients T,-, as we show in Sections III 



IV and [V] For simplicity, we consider an interferometer 



with equal arms, x 1 ^' — x^ 1 = x 1 ^' — x^ 1 = L, which is 
predominantly the experimental situation. In the limit 
E c t 1, where r = L/v is the electron dwell time in the 
MZI, fully analytical treatment is possible. In the more 
general case E c t ~ 1 we have developed a numerically 
exact scheme to evaluate the visibility in the MZI as a 
function of voltage (Se c. |v| ). Before going into details 
of the calculations (Sec. |III| 7we summarize our main re- 
sults. 



(i) 



f>=±' 



which are described by the Grassmann fields ip g in the 
arm q = ±. In the case v > 2 the fermionic fields have 
the vector structure ij) B — (V'lg, • ■ • , ^ug) T due to multiple 
edge channels. The Coulomb interaction is taken into 
account by the electrostatic model with a charging energy 
E c = e 2 /C, such that electrons interact only when they 



B. Limit of strong interaction 

First, we discuss the results in the limit E c t 3> 1. In 
the case of not too low voltages, namely at eVr > 1, our 
model predicts the asymptotic expansion for the differ- 
ential conductance dl/dV of the MZI in the form 




FIG. 2. Scheme of an MZI at filling factor v. Two quantum 
point contacts are characterized by the scattering matrices s 
and s 2 , which connect the outer channels. Inner channels are 
fully reflected. 



+ 2{T 1 R 1 T 2 R 2 ) 1 ' 2 Re 



dip 
d{eVr) 



(6) 



where $ is the magnetic flux and Tq is the amplitude 
of the interference Aharonov-Bohm (AB) contribution to 
the current, 

T = e^+V") (Ci(eyr) Al +G 2 (eV/r) A2 e ±leyT ). 

(7) 

The choice of the sign ± in the exponent of the last term 
will be explained below Eq. (13 1. Equation ^ contains 
the two leading terms of a series (in general, infinite). 
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FIG. 3. Power-law exponents shown as the function of transmission coefficient T\. Solid line shows ReAi, dashed line shows 
Re A2 in the case of v = 3, and (—Re A2) in the case of v = 4, 5. 



The dependence of each term of this series on eVr is 
characterized by a certain power-law exponent A,; and by 
a certain oscillating factor. 

Interpretation of the different ingredients in this ex- 
pression is as follows. The coefficient 



Pi 



1 

2ni 



In Rte 



(8) 



describes the non-equilibrium dephasing of the AB os- 
cillations induced by a combined effect of inelastic e-e 
scattering and the quantum shot noise generated at the 
1st QPC. If v > 3 then Im/3i > and, by defining the 
out-of-equilibrium dephasing rate as 



T7 l = eVImp u 



(9) 



we see that AB oscillations are suppressed by the factor 
e — t/t^ m j-j ic high-bias limit eV S> 1/r. 

It is worth stressing that the exponential suppression 
of the interference current is directly related to the full 
counting statistics (FCS) of electrons passing through the 
QPCf at the time interval r. Indeed, defining the FCS 
cumulant generating function (CGF) of the backscatter- 
ing current as 



Xr(X)= [1 + i?i(e 4A -r 



1 eVr/2-n- 



(10) 



where A is the so-called "counting field'^, we see that 
the damping factor is equal to 



e W = Xr (-Air/v). 



(11) 



The exponents Ai^, which set the power-law depen- 
dence of the interference current on bias, belong to 
the class of non-equilibrium quantum critical exponents. 
Physically, they can be understood as being due to the 
Anderson orthogonality catastrophe which happens each 
time when an electron enters or leaves the interior part 
of the MZI where it strongly interacts with all other elec- 
trons. It is worth mentioning that in the considered sim- 
plified model, where the e-e interaction is present only in- 
side MZI, the orthogonality catastrophe is absent for the 



incoherent contribution to the current, which stays linear 
in voltage as in the case of non-interacting fermions. 

The exponents A 1.2 are functions of both, the filling 
factor v and the transparency T\ of the first QPC, and 
are shown in Fig. [3| The explicit expressions read 



Al,: 



I 1 



2 '^2 



1\ 2 



(12) 



2 < v < 4, 



for low filling factors and 



Ai = 
A 2 = 



2 2 

v v^- 
I 2 

~ 2 + 



(13) 



v > 4 



in the case of higher v. In the case 2 < v < 4 the voltage 
dependent phase factor in Eq. ^ has to be taken with 
the sign (+). For v > 4 the ± sign corresponds to the 
case T\ > 1/2 and T\ < 1/2, respectively. The coeffi- 
cients Ci t 2 in Eq. ^ are some bias-independent complex 
numbers which depend solely on v and T\ and can be 
found from the fit of this asymptotic expansion to its 
numerically exact counterpart. In the limit of strong in- 
teraction, E c t 3> 1, the case v = 1 is very special. Specif- 
ically, one has then I = (eW) and the MZI behavior is 
the same as in the absence of e-e interaction. 

Experimentally, one usually quantifies the coherence 
of the interferometer in terms of the visibility V and the 
phase (Xab of the AB oscillations of the conductance. 
The visibility is defined as the ratio of the amplitude of 
the AB oscillations to the mean value of the conductance. 
In our model 

dl 2(T 1 R 1 T 2 R 2 ) 1/2 , , 

d[eVr) TiR 2 + T 2 Ri 

where Vo is the non-interacting value of V, and 

a AB =a,xg[dl /d(eVT)}. (15) 

In terms of the above quantities the conductance Q(V) 
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FIG. 4. The voltage dependence of the visibility (normalized to its non-interacting value; left) and of the phase (right) of AB 
oscillations for v — 2,3,4, 5 in the strong-interaction limit, E c t 3> 1. Solid lines show numerically exact results, dashed lines 
represent analytical results. The latter are strictly speaking valid in the asymptotic high-voltage limit eVr ^> 1 but turn out 
to work almost perfectly already at very modest values of voltages, eVrf-K > 0.5. 



takes the form 

2 

G(V) = ^(RiT 2 + R 2 T^(l+V(V) cos[<i>+a AB (V)]y 

(16) 

In Fig. [2] we show the visibility V (normalized to its non- 
interacting value Vo) and the phase «ab of the AB os- 
cillations as functions of bias for different filling factors 
v = 2, 3, 4, 5 and for different transmissions Tj. 

In each plot we have fitted the exact visibility (ob- 



tained numerically) by its analytic form based on Eq. 
with two free parameters C\ and G%. Although Eq. |7| is 
strictly speaking an asymptotic formula valid in the high 
voltage limit eVr ^> 1, we see in Fig. [4] that the analyti- 
cal result is an excellent approximation already starting 
from very modest values of voltages, eVr/ir > 0.5. For 
still smaller voltages, the visibility saturates at its non- 
interacting value Vo- 

The most spectacular feature of Fig. [4] are oscillations 



G 



of visibility which become particularly strong yielding a 
"lobe structure with the visibility reaching zero at min- 
ima for v = 2 (for any value of Ti) and for T x = 0.5 (for 
any v). In this cases, the cusps in the visibility at its min- 
ima are accompanied by 7r-jumps in the phase «ab- As 
we explain below, the special role of v = 2 is a character- 
istic feature of the strong-interaction limit. On the other 
hand, the point T\ = 0.5 remains special for a generic 
interaction. 

At v = 2 we have Ai = A 2 = and C\ = -C 2 . This 
gives the oscillatory visibility V = Vo cos(eW/2) which 
is independent of the transparency T\ and does not de- 
cay with bias. The behavior of the MZI in this case is 
analogous to that in a model with short-range electron 
interaction which is also exactly solvable at v — 2 by 
means of the method of refermionization, as has been 
recently shown in Ref. [2SJ On a technical level, the ab- 
sence of dephasing and independence of visibility on T\ 
at v = 2 in the limit E c t ^ 1 comes from the fact that 
the counting phase becomes — Air/v = —2ir in this case. 
At a moderate charging energy E c t ~ 1 both the dephas- 
ing and the dependence on T\ in the visibility of v = 2 
MZI are restored, see Sec. |IIC| 

In the case v = 3 an infinite number of lobes is ob- 
served. As discussed above, the visibility reaches zero 
at minima when (and only when) the transmission is 
T\ = 1/2. The reason for this special role of the point 
T = 1/2 is as follows: in this case the real parts of the 
two exponents are equal, ReAi = ReA2- 

For v > 4 our model predicts only one central and 
one side lobe. Note, that at v = 4 the exponents Ai_2 
logarithmically diverge at T\ — > 1/2. This is the reason 
why at v = 4 we have chosen to plot V(V) and C(ab(V) 
for a slightly different value T\ = 0.45 (see Fig.|4|. 

We turn now to the effect of a dilution of the im- 
pinging current due to the electron scattering at an ad- 



ditional quantum point contact, QPCO, which is put 
outside of the interferometer. At Rq < 1 the QPCO 
generates the double-step distribution function /+(e) = 
Tq9{— e) + Ro8(eV — e) for incoming electrons, which af- 
fects the power-law exponent and serves as an additional 
source of dephasing. The effect of QPCO is particularly 
noticeable in the strong-interaction model with v = 2, 
since in this situation no dephasing and no power-law de- 
cay of oscillations is found in the absence of QPCO (see 
above). In Fig. [5] we show the visibility in this situation, 
with half-transmitting QPC1, T x = 1/2, and for differ- 
ent values of the reflection coefficient Rq of QPCO. In 
the case Rq > 1/2 the suppression of visibility with volt- 
age can be roughly characterized by the dephasing rate 
I/t^ = (eV/2ir) ln(2i? — 1), which diverges logarithmi- 
cally at Rq — > 1/2. At this value of Rq the behavior of 
the MZI visibility changes from the regime with multiple 
side lobes, characterized by periodic oscillations in V(V) 
with a typical period ~ 27t/t, to the regime with only 
one node. Such a transition in the behavior of visibility 
under variation of Rq has been first found in Ref. [32] in 
the weak-tunneling regime 7\ <C 1 for the short-range e-e 
interaction model. Recently, such an effect of QPCO on 
the visibility was observed experimentally^. We believe 
that the experimental conditions of long-range interac- 
tion and Ti w 1/2 are closer to the ones studied within 
our model. 

As discussed in more detail Sec. |IVC~3] the appearance 
of the visibility fringes in our model stems from the su- 
perposition of two multi-particle amplitudes having the 
relative phase shift eVr. These amplitudes correspond 
to processes with different number of particle-hole exci- 
tations between two Fermi edges. In other words, the 
"lobe" pattern in the visibility is a many-body effect 
linked to e-e interaction. 
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FIG. 5. Visibility of AB oscillation at v = 2, li = 0.5, 
and strong interaction, E c t 3> 1, in the presence of an ad- 
ditional quantum point contact, QPCO. The curves from up 
to down were evaluated numerically for the reflection coeffi- 
cient of QPCO equal to R = 0.9, 0.7, 0.5 and 0.3. 



C. The case of moderate strength of interaction 

In this subsection we discuss the results for visibility 
in the case of not too strong e-e interaction, E c t ~ 1. 
We consider here only the case v = 2 for the following 
two reasons. First, the majority of experimental data for 
MZIs has been collected for this filling factor. Second, 
contrary to higher filling factors, at v — 2 a finite (rather 
than infinite) value of E c t changes the result qualita- 
tively, since there is no dephasing and no Ti dependence 
at E c t — > oo. 

The numerically calculated plots of visibility for trans- 
parencies T\ = 0.5 and T\ — 0.2 are shown in Fig. [6] 
It is seen that the finite charging energy E c gives rise 
to the decay of visibility V(V) with bias contrary to its 
behavior at E c t — > oo discussed in the previous subsec- 
tion II B Note also that nodes (zero- visibility points) in 
V(V^) are generally present only in the case T\ = 1/2. 
At the transmission coefficient close (but not equal) to 
one half the nodes are superseded by deep minima. Fur- 
ther, the period of oscillations increases with decreasing 
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FIG. 6. Bias dependence of visibility of MZI with v — 2 for the moderate interaction strength (lu c t = 2, solid blue curve) in 
comparison to the limit of strong interaction (uj c t — > oo, dashed red line), where w c = vE c /-k is the charge relaxation frequency. 
Note that at finite ui c the nodes in visibility are present only at T\ — 1/2. 
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FIG. 7. Time-dependent "counting" phase for v — 2 and two strengths of Coulomb interaction, uj c t = 2 and ui c t = 10, where 



lu c — uE c /n is the charge relaxation frequency. The phase approaches the "window" function, Eq. ( 19 1, in the strong-interaction 
limit E c — > oo. 



E c . However, the estimate e(AV) ~ 2tt/t for the scale 
of oscillations remains valid up to the moderate charging 
energy E c ~ 1/r. As an example, Fig. [6] shows that at 
E c t = 7r the period is larger than its strong-interaction 
limiting value e(AV) = 2tt/t by a factor ~ 1.5. 

The dephasing rate l/r^ describing the exponential 
suppression (oc e _r / r *) of the visibility with bias is found 
to be 



eV 
~2ttt 



rr+t 




/ Re 


In ( 


' — OO 





where t is the time when the electron enters the interfer- 
ometer and 9+ (t) is the time-dependent "counting" phase 
given by 



-(*) 



1 



Im [J> (t- t) - J>{t -t + t)] 



(18) 



with the function J > (t) defined below in Eq. ( |5T| ). The 
phase 0+(t) is shown in Fig. [7] for ui c t — 2 and uj c t = 
10. In the limit E c t 3> 1 the time dependence of d+(t) 
approaches the "window" function 



-(<) 



-n/v, 
0, 



t£ [t,t- 



(19) 



causing the dephasing rate to vanish at v — 1,2. If one 
introduces the effective charge e*(t)/e = 0+(i)/7r, then it 
can be physically interpreted as the optimal charge fluc- 
tuation on the upper arm of the MZI which promotes 
scattering of the transport ("trial") electron from one 
arm of the interferometer into the other Loosely speak- 
ing, if such scattering event starts at a time instant i, 
then it finishes no later than i + r (cf . the upper bound 



of the time integral in Eq. (17)). It means that an elec- 
tron entering the MZI at time t cannot be influenced by 
those electrons which enter at times larger than t + r, 
since by the latter time the trial electron leaves the inte- 
rior interacting region of the system through the second 
QPC. On the other hand, and perhaps somewhat coun- 
terintuitively, a typical arm-to-arm electron scattering is 
generally preceded by a rearrangement of the charge e* (t) 
on the MZI at all times t < i . We thus see that the single 
electron transfer through the MZI in the presence of e-e 
interaction is a collective many-body process involving 
many electrons. 

Our results match well experimental observations in 
many designs of Mach-Zchnder interferometers at fill- 
ing factor v—2, which happen to be rather universal. 
Namely, at T\ close to 1/2 the experimentally observed 
dependence of the visibility on voltage shows a number 



of "lobes" whose amplitude gets suppressed with the in- 
crease of bias. At the same time, the voltage dependence 
of the AB phase is close to a piecewise constant func- 
tion with jumps of a magnitude 7r at minima of the vis- 
ibility. Further, we estimate the period of oscillations. 
As follows from Fig. [6j the characteristic energy scale 
corresponding to the first minimum in the visibility is 
e(AV) ~ 2it/t — 2irv/L. An estimate for the drift veloc- 
ity in our phenomenological model, v ~ ve 2 /eirh, can be 
obtained following Ref.|54j where the excitation spectrum 
of the compressible Hall liquid has been studied (see Sec. 
IV. C of our previous work, Ref. |351 for a more detailed 
discussion). Taking e = 12.5 for the dielectric constant of 
the GaAs heterostructure, one obtains v ~ 1.1 • 10 5 m/s. 
Note, that this estimate agrees well with an effective ve- 
locity v c s = 6.5 • 10 4 m/s found in Ref. HOI from the anal- 
ysis of data on energy relaxation in QH edge states at 
v=1. For a typical size of the interferometer L ~ 10 /im 
we then get AV^ ~ 40/iV, which is of the same order as 
the experimentally observed energy scale of the visibility 
oscillations. 

Having completed a presentation and discussion of our 
key results, we now turn to the exposition of the method 
and of technical aspects of the derivation. 



III. NON-EQUILIBRIUM FUNCTIONAL 
BOSONIZATION FOR MACH-ZEHNDER 
INTERFEROMETER 

In this Section we show how the method of the non- 
equilibrium functional bosonization can be used to solve 



the model of the MZI defined in Sec. II A First, we 



present the Keldysh action of the problem and derive 
the expression for the direct and interference current 
(Sec. |III A ). Then we give details of the real-time non- 
equilibrium instanton approach (Sec. IIIB). Using the 
special structure of the Keldysh action, we show that this 
method becomes exact in the case of the simplified model 
of the Coulomb interaction (considered in the present pa- 
per) in which electrons interact only in the interior region 
of the MZI. Finally, we specify the form of the instanton 
for the case of the constant interaction model. 



A. Keldysh action and current 

The theoretical model of the MZI, which we consider 
throughout the paper, is defined by Eqs. and To 
make our discussion more general, we will first assume 
the arbitrary interaction potential U± (x — x') between 
two electrons in the same edge, which however is non-zero 
only if 2 ]. Because of the non-equilibrium 

character of the prob lem, we proceed within the Keldysh- 
type frameworlPHH!. We decouple the interaction term 
iSjnt using the Hubbard-Stratonovich transformation with 
fields f e (x, t), where the index g = ± labels two arms of 
the interferometer. Following the logic of the Keldysh 
formalism, we then double the number of Grassmann 
fields, ip e = (ipg, ipg), as well as of the bosonic fields 
ip B = (ipfcp), where indices / and b denote the fields 
residing on the forward and backward branches of the 
Keldysh contour C, respectively. These steps lead us to 
the MZI action in the form 



S = / dtdxip e (idt + ivd x — tpg)ip e + - / dtdxdx' (p g (x)U e 1 (x — x')<p e (x'). 
e =± J c - h 



Integration along C is to be understood as J c dt'A(t') = 
fdt>(Af(t>)-A»(t>)). 

In terms of fermion fields ip g the action S is quadratic, 
thus they can be integrated out. In this way we obtain 
the Keldysh action A[ip] of the MZI which depends on the 
electrostatic potentials <p±(x,t) on two arms. The out- 
lined method is known as the functional bosonization. 
The integration over the Grassmann fields ipg should to 
be performed with taking into account the relation Q 
at QPCs; this relation has to be satisfied independently 
on each branch of the Keldysh contour. The action A[<p] 
in the case of a generic non-equilibrium setup formed by 



ID electronic channels coupled by a number of local scat- 
ters and by electron-electron interaction ( "quantum wire 
network") has been found in our previous worlP^ In 
particular, the Keldysh action A[(p] describing the MZI 
can be expressed in terms of the time-dependent single- 
particle scattering matrix of the interferometer in the 
given configuration of the fields ifi^ b , which we denote 
as Sf/b — S[tpf/ b ](t, t'). This S'-matrix describes electron 
scattering at both QPCs and the propagation of elec- 
trons along the arms of the MZI. The bosonized action 
A[4>] has the form A = A, n t + -4f cr m + SA. Explicitly, 



Ant = 2 



E 

Q 



« <P c e (Z) {UgHx, x')5(t - o - n£(£, a) *%(£), 



(20) 
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where we have denoted £ = (x,t), and 

A form + <5 A = —iTi In 



1-f + Sle^Sff] -2^Tn?g e / ( 



(21) 



Let us explain notations and comment on different terms 
in the action, Eqs. (20) and (21 ). In the above expression 
for Ai n t we have used the Keldysh basis ip g = (ip c g ,ip q g ), 
with "classical" and "quantum" components being de- 



fined as ipg ,q 



(<pf ± <p h g)/2- The advanced component 
H A of the ID polarization operator, when written in 
the frequency-coordinate representation, has the explicit 
form 



U a (uj,x,x') = 



2ttv 



S(x - x') (22) 



l-e-'*" 1 ''/' x 9{x-x') 



The combination of the polarization operator and the 
bare interaction potential entering Eq. ( 20 ) determines 
the RPA screened interaction potential, 



v: 



-1,A 



U- X (x,x')5{t- 



t') 



n^,0- (23) 



The action .4f orm has the form of the fermion deter- 
minant and bears a close connection with the problem 
of electron full counting statistics^. We have intro- 
duced the distribution functions of the source reservoirs, 
/ = diag(/ + , /_), where f± are diagonal matrices in the 
channel space. Without QPCO, To = their components 



are/5(t,t / ) = e 



_ _-ieV£(t-f) 



fo(t-tf). Here V£ is the volt- 



age applied to the n-th channel in the upper/lower arm 
and 



fo(t-t') 



1 



2nt-f 



iO 



(24) 



is the time representation of the equilibrium Fermi distri- 
bution function. For instance, for the MZI presented in 
Fig. [l]in the case of vanishing transparency T = 0, the 
only non-zero voltage is Vl = V. At To > 0, when QPCO 
is used to dilute the impinging current, the function /j. 
is the double-step distribution in the energy domain and 
is given by 

f\{t - t') = [j2o e~ ieV ^ + T ] f Q (t - t>) (25) 

in the time representation. We have also introduced 
the auxiliary "counting fields" \ — diag(x + ,x_) in 
the drains which enable us to find the number of elec- 
trons transferred through the MZI. The determinant in 
Eq. ( |2l| ) is evaluated with respect to time, channel and 
arm indices. 

Next, we specify the S'-matrices Sf/f, of the MZI which 
enter Eq. ( 20 ) and encode all information about electron 



scattering. We introduce the phases d^J h (t) accumulated 
between the QPCs 1 and 2 due to interaction, 

H /b (*) = / dx' d /b (x\ t+(x'- xl)/v). (26) 



Let now x g and x® be the coordinates of the source and 
drain reservoirs. We also define the "time-delay" opera- 
tor A lk = diag(A+, A^), where 



A l e k (t,t') = 5(t-t'-(x l e -x k e )/v), 



(27) 



with indices k, I G {S, 1, 2, D}. It coincides with a trans- 
fer matrix from x k e to x l g along the arm g of the MZI in 
the non- interacting limit. Assuming the absence of inter- 
action outside of the interferometer, we obtain the total 
5-matrices 



5' 



f/b 



A D2 S 2 A 



21 ?:* ii9 



Si A 



is 



(28) 



Here §i and S 2 are the local scattering matrices of the 
1st and 2nd QPC. Further, <J> is the diagonal flux ma- 
trix $ = diag(<I>/2, -$/2) ffi t2u-2, where $ denotes the 
Aharonov-Bohm phase and the sign ± distinguishes be- 
tween the upper and lower arms. The direct sum (©) 
refers to the channel space and l n is the n x n- unity 
matrix. The matrix df/ b = diag(?9 



f/b fif/b 



+ 



has an anal- 



ogous structure. For the MZI scheme shown in Fig. [2] 
only the outer channels are mixed by scattering. In this 



case, Sj — s J © t 2u -2, with s J given by Eq. rt5 



Finally, the counter-term 8 A in the action (21) is in 



eluded to cancel the equilibrium Fermi-sea contribution 
which does not affect the non-equilibrium electron trans- 
port. The "quantum" phase i?Q g entering S A is defined 
as 0g e (t) = df(t + (xf - x\)/v) - d e {t + (xf - x\)/v)- 
the trace (Tr) is taken over channel and arm indices and 
also includes integration over time. 

The bosonized action A[(p, x] enables us to find the 
generating function of the interferometer's FCS as a func- 
tional integral over ip, 



Z(x) = / V^J\x,t) wp{iA{<p,$)}. 



(29) 



Then the number of electrons transferred to, say, the 
lower drain during a long observation time to 3> 
max{h/eV,hv/L} is obtained as 



N- = -id x _ InZ = ( d x _Af etm 

x=o \ x=0 



(30) 



The quantum mechanical average (. . .) v here is under- 
stood as the path integral over <p with the weight e lA , see 



Eq. ( 29 1 , but with "counting fields" \ P u * to zero. Since 



the Coulomb interaction is assumed to be absent outside 
the interferometer cell, A{ e rm simplifies considerably (in 
the rest we will not explicitly state \ — anv longer): 



A {e 



-ilnDetD, V = t 



f + S\e 2i ^ St f. (31) 
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The same action can be represented in the equivalent 
form as A[ CTm = — ilnDetP, with 



V = S!VSl = t-f + e 2M f, 



(32) 



where / plays a role of the non-equilibrium density ma- 
trix of the interferometer. If the voltage is applied to the 
outer channels only, then / — f 1 (I21/-2 • Jo), with 



I 1 = hM = 



-iy/R^Ul-f-) 



(33) 

Note, that -4f crm at \ — depends on the scattering 
matrix s\ only. It also depends solely on the "quantum" 
field ip q , but not on the "classical" one. These special 
features stem from the chiral nature of the MZI. The 
independence of *4f crm on the classical component of the 
field will play a crucial role in the sequel, as it will allow 



us to find an exact solution of the problem. 



Using the definition (30) and the full ^-dependent 



fermion action (21), one obtains the following interme- 



diate expression for N- 



.V_ =('&X>- 1 ate-*^ e - < *A 2lt 4| 



(34) 



To derive this result we have taken into account that 
only outer channels with a non-equilibrium distribution 
matrix f 1 may contribute to the transport of electrons 
and have introduced the basis |±) in this linear subspace. 
Taking the trace over the c han nels and using the explicit 
form of A 12 , given by Eq. (27), we find N_ — J2 
where 



N 

flK 



N, 



(IK 



! s 2 \iA 



K\ SA 



dt 



(l4 «1 (/ 



o2T 



(35) 



In this expression we have introduced <E>± = ±$/2, which 
are just the AB phases accumulated at each arm of the 
MZI, and defined t± = (x\ — x±)/v — the flight time of 
electrons along upper/lower arm. The sign "o" denotes 
the convolution in the time and channel space. Clearly, 
the diagonal (N^) and off-diagonal (iV^-^) elements 
give, respectively, the direct and interference contribu- 
tions to the total current. 



B. Exact solution: from many-body to 
single-particle problem 

In general, the functional integral for interacting elec- 
trons in a quantum wire network cannot be evaluated 
exactly. In Refs. [35l [50J, and [57] an instanton approach 
has been developed which yields a controllable approx- 
imation to the problem for the case of weak tunneling 
between the channels. It turns out that for the problem 
considered here this method becomes exact (for any tun- 
neling strength). Specifically, we show below that the 
functional integral can be exactly evaluated in a fash- 
ion similar to the exact solution of the problem of non- 
equilibrium Luttinger liquids in Refs. HH [45j and EH In 
fact, we will see later that there is a deep connection 
between the two problems. 

We have shown above that the number of electrons 



expression implies the path integral over all realization of 
the fields ^\x,t). As we reveal below this functional 
integral can be performed exactly. What crucially sim- 
plifies the calculation of the quantum mechanical average 
is the fact that T> and hence the /^K-matrix elements in 
Eq. ( 35 ) together with Af clm do not contain ip c or, equiv- 
The classical fields enter the RPA action ^4j nt 
and appear there only 



alently, $ c 



and the phases $ f e /b = & c g ± 
linearly. Therefore ip c can be exactly integrated over. To 
this end let us introduce sources J so that 



A 



J;fiK 



^(i)-^(i + T„-T«) 



(36) 



We see that for // = k the source vanishes. On the con- 
trary, at \x = — k the explicit expression for J reads 

J q ^ K (t,0 = ~v- 1 S(t + (x - xl)/v - t), (37) 
J 9 ^ K (t,0 = v- 1 5(i+ t^-t k + (x - xl)/v - t). 

The physical meaning of the source terms in the action 
is rather obvious. They describe an electron transfer be- 
tween two Hall edges of the MZI, thereby creating a hole 
in the arm k and adding an e xtra electron into the arm /i. 
Decomposing i}^ b in Eq. (35) into the "classical" and 



transferred through the MZI is given by Eq. ( 35 1 . This 



"quantum" parts, we rewrite the formula for the particle 



numbers A^„ K in the form 



dt 



JAf, 



-,+iSA 



Afj,K (t>t 



r«)}, 



(38) 
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where the prefactor Auk is defined as 



(39) 



As has been emphasized previously, the fermion action Ai etm + SA and the matrix T> are Junctionals of tp q only, 
see Eq. (31 1. Hence, one can first perform the integration over the "classical" field ip c . Taking into account that 



is linear in <p c , we obtain 

J V<p c exp [iA int + iAj.pt] = J V V C exp [2i V c {V A )~ V + ^ c J q \ <x6(<p«- (pi) . 
where the (5-function fixes the quantum component tp q to be equal to the saddle-point trajectory 



(40) 



(41) 



The (5-function constraint renders trivial the subsequent integration over (p q . Taking quantum-mechanical average 
(. . .) v is therefore reduced to the evaluation of the integrand e lAtt2TTCl+lSA A^ K on the optimal trajectory ip q = ipl. The 
particle numbers are then simplified to 



N 11k = (-\s 2 \h)(k\sI\-) / die 



(42) 



with the "quantum" saddle-point phase (or "instanton" ) 



A^ =Rifl + Txf^, that yields 



n e {t) = -v- 1 / dx' V Ux', t+(x>- xl)/v e ). (43) 



In what follows we will frequently refer to M as to "count- 
ing phase" in view of an analogy between the action Ai erm 
and the theory of the FCS. 

To reiterate the logic, we have reduced the path inte- 
gration over tp to the evaluation of the integrand for the 
numbers N„ „ on the "quantum" saddle-point trajec- 
tory ip%. This is the main result of the present subsec- 
tion. The optimal "quantum" electrostatic field is re- 
lated via the RPA interaction potential to the source 
terms J which describe the electron transfer between 
two edges of the MZI, see Eq. pi] ). The bare interac- 
tion potential enters the result through the RPA kernel 
V" A (£,£'), thus the outlined method is very general, pro- 
vided the e-e interaction can be disregarded outside of 
the MZI cell. The result is expressed in terms of deter- 
minants and resolvents that are of single-particle com- 
plexity; thus, we have achieved a dramatic simplifica- 
tion as compared to the original many-body problem. 
Needless to say, the obtained apparently single-particle 
quantities carry all the physical information about the 
many-body physics of the problem, including, in par- 
ticular, non-equilibrium orthogonality-catastrophe expo- 
nents and non-equilibrium dephasing. 

Let us now discuss the direct (incoherent) contribution 
to the current which arises from the diagonal numbers 
{N ++ and TV ). Within our model, they are not af- 
fected by e-e interaction. Indeed, in this case the sources 
vanish, J^. MAl = 0. Hence the instanton trajectory is triv- 
ial, §1 = 0. One thus get T> = 1, „4f orm + 8 A = and 



N++ = T 2 Tr [ J R 1 /++T 1 /_] 
N__=R 2 Tt [R x f_+T x f+] 

Taking the Tr in the time space, we have 

eVt 



(44) 



N++ + JV_ 



2nh 



(l^x+^Tx). (45) 



Thus, the direct current is linear in bias; it is the same 
as in the non-interacting limit. 



C. Constant interaction model 

In this subsection we specify the "counting phase" for 
the constant interaction model with the charging energy 
E c , as it is defined by Eq. ([2]). We assume that both arms 
of the MZI have the same length L, thus t = t + = t_ 
and x 1 ^ 2 = x^ 2 . Then the bare interaction potential 
U±(x, x') = E c if both x, x' £ [x 1 , x 2 } and U± (x, x') = 
otherwise. This form leads to the RPA potential V A {ui) 
which is non-zero only if X . X ^z. ^X - X 1 and constant 
inside this region, 

V^(u))=E c (l-E c J dxdx'n A {uj;x,x')j . (46) 



Using Eq. ( 22 ) one has 



dxdx' n (w; x, x') — 



iv 1 — e 

277 U) 



(47) 
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Defining the charge relaxation frequency as ui c 
vE c /(2tt) we obtain 

V A {u>) = ^ 



iuJl 



1 



(48) 



We now use this RPA result to find the instanton po- 
tential ipl and the phase d%. By virtue of Eq. (411 we 
have 



2T 



dx' V A {t-t- (x'-x^/v), 
(49) 

where lr x i )X 2i(a;) is the projector on the interval where e- 
e interaction and thus the potential are present. Taking 
into account the relation ( 43 ) between the phase and the 
potential, one finds 



dx'dx" 



du> 



= ^Im \j>(i 
v L 



2tt 

t) - J>(T-t 



"y^V A (u) 
t)] , (50) 



where we have introduced the phase-phase correlation 
function 



J>(t) 



did 



u(uj + iu} c (l - e iulT )) 



1) 



The function J > (t) has appeared and was analyzed in our 
previous work (see Sec. 4.3.4 of Ref. IST)|) in the context 
of theory of Fabry-Perot QH interferometer. In the limit 
of strong coupling u> c t 3> 1 and long time u c t 3> 1 it can 
be well approximated by the logarithmic asymptotic 



J>(t) 



"7 



In 



t + ia 



(52) 



with a being the short-time cut-off, sec Fig. U Therefore, 
except for times f in a close vicinity of either t or t + r, 
the "counting phase" simplifies to 



(53) 



with the K-dependent constant i?^ and the unit "window" 
function 



mAt) = d(t-t)-6(t + T- 1). 



(54) 



In the case of a moderate charging energy, uj c t ~ 1, one 
has to resort to a numerical evaluation of the (imaginary 
part of) correlation function J > (t). We have — 
—Krj6+(t), where the function 9+{t) is independent of k 
and rj. Typical plots of 6+(t) (at v = 2) are shown in 
Fig. [7] of Sec. |II C| For brevity we will omit the * index 
when denoting the counting phase ^{t) = J?*^(t) in the 
following. 

In passing we note, that the time-dependent counting 
phase i}+(t) in our theory is the analog of the kernel Q(x), 
introduced in Ref. [33] (cf. Fig. 7 in this work). Similar 
to the phase #+(t), this kernel depends on the nature of 
interaction potential and is used to describe the phase, 
which an electron, passing the MZI, accumulates due to 
interaction with other electrons. 
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FIG. 8. Correlation function J > : comparison of analytic ap- 
proximation J^ a given by Eq. ( 52 1 and numerical results J, 
for ui c t — 25. 
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STRONG INTERACTION: ANALYTICAL 
SOLUTION 



IV. 



As shown in Sec. |III| the considered model of a MZI 
with inside-only interaction can be exactly solved by the 
non-equilibrium functional bosonization method. The re- 
sult is expressed in terms of single-particle objects: de- 
terminants and resolvents of Fredholm operators. While 
it is not too difficult to evaluate such quantities numer- 
ically, it would be highly advantageous to have a fully 
analytical solution of the problem. Such a solution will 
be obtained in the present section for the regime of strong 
interaction, E c t 3> 1. This solution will allow us to un- 
derstand much better the physics of the problem, includ- 
ing the formation of the visibility oscillations (taking a 
form of "lobes" in certain situations as discussed above) 
and their decay with voltage. Further, while determining 
the exact solution, we will establish a deep connection of 
the present problem with that of non-equilibrium Lut- 
tinger liquid and, more generally, with a broader class of 
non-equilibrium many-body problems. 

As we have demonstrated in Sec. 
phase" in the strong-interaction regime 
duced to the "window" function on the interval [t, i + r] . 
We will show in Sec. IIV Al that under this condition the 
interference current can be represented in terms of sin- 
gular Fredholm determinants generalizing Toeplitz deter- 
minants with Fisher-Hartwig singularities. Using asymp- 
totic properties of such determinants (Sec. IV B I, we fur- 



III 



the "counting 
E c t ^> 1, is re- 



ther derive the high-voltage form of the AB contribution 
to the current (Sec. |IVC]. The result is Eq. (171 which 



has been already discussed in detail in Sec. |IIB 



A. Reduction to a single-channel problem 



In Sec. IIIB we have expressed the interference current 
in terms of the operator T> and the non-equilibrium den- 
sity matrix / which, in addition to being the matrices in 
the time space, have also a nontrivial channel structure: 
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for given times t\,t% they are matrices from C 2x2 . (Since 
all the relevant non-equilibrium physics arising due to 
scattering at QPCs concerns only the outer channels, we 
focus on these two channels. In what follows we con- 
sider projections of all operators, such as /, onto the 
two outer channels; thus they retain the smaller 2x2- 
channel structure.) The double index structure (times 
and channels) very seriously complicates the computa- 
tion of the determinant and finding the inverse of T>. In 
this section, using the Riemann-Hilbert technique, we re- 
duce the matrix determinant and resolvent to a product 
of certain single-channel determinants. We show that 
the corresponding operators belong to the class of sin- 
gular Fredholm operators that may be considered as a 
generalization of Toeplitz matrices with Fisher-Hartwig 
singularities. The reduction to a single channel problem 
will allow us to calculate analytically the current in MZI, 
see Sec. MM and MH below. 



1. Heuristic argument: Relation to full counting statistics 

Before presenting a rigorous derivation of the reduction 
formula, we will put forward a more heuristic argument 
in its favor. This argument is based on a connection be- 
tween the fermion action -Afcrm and the theory of the full 
counting statistics (FCS). Consider the cumulant gener- 
ating function (CGF) for the statistics of N±, the num- 
bers of non-interacting electrons which tunnel through 
the QPC1 during the time interval [i, t + t] into the up- 
per/lower arm, resp. 



Xr (A + ,A_) = ( e ^+ + ^-) 



(55) 



with corresponding "counting fields" A±. The brackets 
here mean a quantum-statistical average. It is known^ 
that this CGF can be represented as a functional deter- 
minant, 



Xr(A+,A_) =Det 



1 



f + s\e iX sJ 



(56) 



where / = diag(/ + , /_) is the incoming distribution ma- 
trix of the 1st (outer) channel of the MZI. The counting 
fields here are assumed to have a time-dependence given 
by the "window" function (54 1. In this way the measur- 



ing time is encoded in the above formula. Let /A is set to 
the equilibrium Fermi distribution, and is the Fermi 
distribution with the chemical potential eV. Asymptot- 
ically, at eVr ^> 1, and dropping the equilibrium contri- 
bution i(A + +A_)(t/27t) J* de (which is infinite because 
of the chirality), we obtain 



In kr(A + ,A_)] 



eVr 



ln[i? ie lA + +Tie lA -] . (57) 



By comparing the Eqs. (56) and (31) we conclude that 
in the limit E c t 3> 1 



Mfe 



DetP = X r (2i?+,20_) . 



Because of this relation to the theory of FCS we will 
frequently refer to the matrix T> as the "counting oper- 
ator". In the presence of scattering at QPC1 this op- 
erator possesses a 2 x 2-matrix structure in the channel 
space. Consider further a single chiral channel with some 
(in general, non-equilibrium) distribution function / and 
the phase S(t) = Wt T (t)S. Then the generating function 

A T [<5, /] = (e lSN ) of the number N of electrons passing by 
some observation point during the time interval [t, i+r] is 
given by the determinant of a "scalar" counting operator 
of the kind 

A T [<$,/] = Det [1 + (e l5 - 1)/] . (59) 

We now argue that the determinant of the matrix count- 
ing operator 2?, evaluated on the Fermi-like distribution 
functions (which is the case of reflectionless QPCO with 
Ro = 1) can be factorized into a product of "scalar" de- 
terminants of the type (59). 



Let us take a closer look on the random numbers N± 
in the CGF given by Eq. (55). In the strongly non- 
equilibrium situation which we consider, i.e. when volt- 
age dominates over the temperature, they should be 
significantly negatively correlated (due to partition at 
QPC1), while their sum, N = N + + N- should be much 
less sensitive to scattering and will only weakly fluctu- 
ate around (N) = eVr/(2n) + N with N being some 
equilibrium contribution. (At the strictly zero tempera- 
ture and long time limit N does not fluctuate at all.) We 
thus expect that N and N+ are only weakly correlated. 
Taking further into account that ?9 + = — we obtain 

X r (2^+,2^_) = ( e 2«n* + + 2 "^-) = 



i e 4i-»\N + 



He- 2 " 5 '*). (60) 

This representation maps our 2-channcl problem to three 
single-channel ones. The term (e 4i-l3 + +) requires us to 
count the charge in a single channel of the upper arm. It 
is characterized by the distribution function 

f ++ (e)=R 1 0(eV-e)+T 1 6(-e) 

established by the QPC1, thus we conclude that 



4«9 + Ar+N 



A T [4tf + ,/ ++ ] 



(61) 



The second term, (e" 2 ^+*) = ( e -™+(N++N-)) y counts 
the total charge in both arms. In view of charge con- 
servation this total charge can be also measured in the 
incoming channels the contribution of which are uncor- 
rected. The corresponding distribution functions are 
/+(e) = 6(eV - e) and /_(e) = 0(-e). Therefore, we 
get 

(e-**+*) = A T [-2tf + ,/+]A T [-2tf + ,/_]. (62) 

Combining the relations ( 60|62 ) we arrive at the following 
result: 



X x(2^+,-2^ + ) 



(58) 



= A T [4^+,/ ++ ] 

xA r [-2tf +1 / + ]A T [-2tf +) /4 (63) 
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Finally, we complete the argument by using Eq. ( 58 1 , 



which yields the desired decomposition of the matrix de- 
terminant DctI? into a product of scalar determinants: 



DetV = A r [4tf +,/++] 

x A r [-2tf+,/ + ]A r [ 



(64) 



The following remark is of order here. The equal- 
ities between the determinants arising in our context 
and in the context of the counting statistics are, strictly 
speaking, valid for not too large values of the counting 
phase. For larger phases the counting statistics determi- 
nants show singularities and switch to another Rieman- 
nian sheet, while our determinants behave analytically, 
see Refs. 135] and [55] for an extended discussion. Physi- 
cally, this is due to the fact that the counting statistics 
"knows" about the charge quantization, whereas for our 
problem the charge quantization is of no relevance. This 
remark does not affect the validity of the final result (64 1, 



since both sides of this equation are analytic functions of 
the phases. 



2. Rigorous proof: Riemann-Hilbert method 



We are now going to prove Eq. ( 64 1 rigorously by an- 



alyzing the associated Riemann-Hilbert (RH) problem. 
We consider the function 



Y(t) = exp 



7T t-t- 



(65) 



which is analytic and non-zero in the complex plane t £ 
C, except for the interval of real times [i, t +r]. It has 
the property Y(t) — > 1 at \t\ — > oo. Next we define 
the functions Y± (t) — Y(t ± iO) on the real axis, t £ 
K. They solve the (scalar) RH problem YI 1 (t)Y + (t) = 

e - 2 *»+M, where i9«_(t) = j} + w {>T (t). The functions Y±(t) 
obey important identities, 



/< r_ f< = ft y_ , f> r_ f> = r_ f> , 

/o Y+ ft = Y + f , / Y + f = f Y + , 



(66) 



where the convolution in time on the left-hand side of 
these relations is implied and we set ft(t) = fo(t) and 
fo(t) = &(f) ~ fo(t ) — fo(—t). As an example, the first 
relation in Eq. ( |66[ ) reads in explicit notations as follows: 

i/2w 
h-t 2 + iO 



i/2n i/2n 
h-t + iO 't-t 2 +iQ 



Y-(h). 
(67) 

Due to analytical properties of y_ (t) this integral is de- 
fined by the residue at t = t<z — iO in the lower half-plane. 
We also note that in the energy domain at zero temper- 
ature, fo~ is the projector on occupied states, whereas 
/q projects on unoccupied states. Therefore, we have 

(/(f) 2 = ft, (fo? = fo, fofo = = ft ft, and 
f$ + ft = 1. The same relations hold in the time do- 
main as well, where the product of two operators is un- 
derstood in the sense of convolution. Using the basic 



identities (66) one derives another two useful relations, 

/o < ^-/o > =0, / > r+/o < =0. (68) 

Let us now turn to the analysis of the counting opera- 
tor V defined by Eq. (31). We introduce the gauge ma- 
trix A = diag(e _4el/ +*, e _leV ^*) comprising voltages V± 
applied to the outer channel in the upper/lower arms. 
By using these gauge factors one rewrites T> as 



1 t 2i!9 9 



V = Alf>+A-'s[e 



siA/ < 



(69) 



With the use of solution to the RH problem we have the 
identity (#1 = -01) 



a- 1 s \ 







siA 



A 



1 



s x A Yl x Y+ 



(70) 



Bearing in mind that Y_ 1 is a local in time operator 
without matrix structure in the channel space, one can 



commute it to the left of Eq. (|70| . In this way we find 
V = KYI 1 



Y_f>+K- l s 



1 



siAY+f< 



A" 



(71) 



To proceed further we apply the unitary transformation, 
T> = siT)s\, and factorize the operator T> into a product 
of scalar counting operators. This is possible by virtue 
of the identity 

(ft + Aft) (Y.ft + Y + f<) = Y_ft + AY + ft, (72) 

which is valid for a local in time 2x2 matrix A(t). As 
one can check, the above relation follows directly from 



the projector properties, given by Eqs. (66 1 and (|68j). By 
setting 



.4 



A" 1 s\ 



1 



siA, 



(73) 



we obtain 
V = 3C 1 SiA 



/o > +A"M 



e 4»n o 

1 



siA/ < 



A 



-M 



x Sl A[r_/ > + r + / <] a- 1 ^. 

If one further introduces operators 



v 



+ 




1 




1 /, 



(74) 

(75) 
(76) 



where 



/ = a 1 fa\ 



t _ 



/++ /+- 
/-+ /— 
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is the non-equilibrium density matrix in the MZI cell. 

Explicitly, we have /++ = R\f\ + Ti/i and /_| = 

i(i?iTi)V2(/i _ /i). Then V is equivalently rewritten 



V = Yl l V* Sl Y-V sl. 



(77) 



To obtain the operator T>q we have used here once again 
the solution of the RH problem. We hence conclude that 



Det V = Det 2X Det X>, 



o- 



(78) 



It is now straightforward to evaluate two determinants 
appearing on the right-hand side of this relation. In the 
case of matrix 2X we obtain 



2X 



(79) 



V 1 / 

where the scalar (in the channel space) counting operator 

= I + (e m + - l) /++ (80) 

is expressed solely in terms of the upper diagonal block of 
the density matrix /++, which, obviously, has the mean- 
ing of the non-equilibrium distribution function in the 
upper arm of the MZI. As the result, 



Det £>*=A T [4tf +,/++]. 



(81) 



In the case of matrix T>q the incoming density matrix / 
is diagonal in the channel basis, that yields 

Det V = A T [-20+, /+] A T [-20+ , /_] . (82) 



Combining together Eqs. (78 1, (81 1, and (|82j), we obtain 
the relation ( 64 ) . The proof of this formula is thus com- 
pleted. 



3. Inversion of the matrix counting operator 

In the preceding subsection we have proven that the 
determinant Det 2? (or, equivalently, the fermion action 
-4ferm) can be expressed in terms of determinants of 
single-channel (scalar) operators. For the evaluation of 
the interference current one also needs to consider off- 
diagonal matrix elements (/i = — k) (/x|/2? _1 (t,i)|/c), see 
Eq. ( 39 ). The goal of this subsection is to show that, sim- 



ilar to the action, the above matrix elements can be also 
expressed via the scalar counting operator 2?**, given by 
Eq. pOl. 



According to Eq. ( 77 1 , the inverse of V can be written 



as 



t)- 1 = s 1 Vq 1 YZ 1 s{v- 1 Y_ 



(83) 



Making use of the solution to RH problem we then rep- 



resent the counting operator ( 76 1 in the form 



V Q = AFT 1 (F_/> + Y + f<) A" 1 . 



(84) 



The basic relation of the RH method, 

(y_/ > + Y + f<y' = Y = \f> + Y^f< 



(85) 



easily gives the inverse of 2? (one can check the former 
identity by multiplying two operators to get the unity, 
employing for that relations (661 and (|68|). The inverse 
of T> then reads 

V- 1 = Sl A [Y- X f> + Y^f<] k-hlv-'Y., (86) 

and the subsequent convolution with the MZI's density 
matrix / yields 



fV- 1 = Sl A/ A 



Y^fV-'Y. 



(87) 



The required (/it, nt) matrix element of this operator then 
takes the form 



e 2 ^ M/B-^t-t)!*). 



The inversion of the operator appearing here is not 
exactly trivial, but it is simplified a lot due to its tri- 
angular structure (79) in the channel space. Note that 
the relation X>**(ti,l2j = 8(ti — t 2 ) for t\ ^ [i,t + t] im- 
plies the same for the inverse, V~^{t\, t 2 ) — S(ti — t^) for 
t\ £ + t] (this can be seen by employing the block 
matrix representation or using the reformulation in terms 
of the Riemann-Hilbcrt problem) . One therefore obtains 

(-| fV-\t,f) |+) = /_+2?- 1 (t,t) = f-+w lT V-}{lt), 

(89) 

and the analogous relation for the conjugated matrix el- 
ement, 



= DJwlt f- 



-(*,*) 



1 /+- 



t,t 
(90) 



It is worth pointing out that the instanton phases = 
ztrjir/iy in Eqs. (89) and (90) have opposite signs (and 
hence 2?* and 2?** differ between these two equations). 
Since under complex conjugation f v (t)* = f v (— t), these 
two matrix elements are indeed complex conjugates of 
each other. Relations (89) and (90) are the final result 



of this subsection and will be used below in Sees. II VI and 
|V]foi" evaluation of the interference current. 



B. Toeplitz matrices and their generalizations 

In this subsection we relate the current and the action 
to the theory of Toeplitz matrices. We review key results 
on the large- N asymptotic behavior of Toeplitz determi- 
nants with "Fisher-Hartwig singularities" and of a more 
general class of singular Fredholm determinants. These 
results will be then used to calculate the determinant and 
the inverse of the operator 2?**, which will serve as the 
basis for the calculation of the AB conductance made in 
Sec. ESI 
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1. From integral operators to Toeplitz matrices 

We relate first the fermion action Af erm to the theory 
of Toeplitz matrices. As we have shown in the previous 
section, the action is expressed in terms of single-particle 
counting operators, see Eq. (64). 



The linearized single-particle spectrum used in our ID 
model lacks upper and lower band edges. Thus, a defi- 
nition of the determinant of such operators requires an 
ultra-violet (UV) regularization. One possible way to im- 
plement the regularization is a discretization of the time 
coordinate t in steps At = 7r/A, which amounts to the in- 
troduction of an UV cutoff A and restriction of energies 
to the range [—A, A]. In this regularization procedure 
operators with kernels such as X>** (£i , £2), cf- (80), are 



turned into (in general, infinite) matrices with discrete 
time indices. 

In the limit of strong interaction the phase = d v Wt tT 
is a piecewise constant function which vanishes outside 
the interval [t, i + r] . Introducing the projector P which 
acts on functions (f>(t) by multiplication with a window 
function in time, P<j> = <j>Wt T , and thus satisfies P 2 = P, 
we can write 



Det X>* 



Det 



t + P 



1/4 



= Det 1 + P ( e 4l1? + - 1 ) J++P 



(91) 



The operator gr„ = 1 + P (e il ®+ - l) f ++ P has a block 
structure, namely 



g**(ti, t 2 ) 



_ / g(h -£ 2 ), h,t 2 e [t,t + r] 

<5(£i — £2), otherwise, 



(92) 



where 



g(ti - h) = S(h - t 2 ) + (e 4W ~+ - l) /++(£! - t 2 ). (93) 

The kernel of the operator is nontrivial only if both t\ 
and t 2 lie in the interval [£, t+r], in which case it depends 
solely on the difference t\ — t 2 . The determinant of g** 
will be given by the nontrivial block g{t\ — 1 2 ). The UV 
regularization of g as described above will give rise to a 
large A x A-matrix (ftfe)i<i,fc<An A = r/Ai = At/tt, 
whose elements depend on index differences, gjk — gj-k- 
The matrix (gjk) of such type is known as a Toeplitz 
matrix. 

Matrices of Toeplitz form are ubiquitous in mathemat- 
ics and physics where they appear in a variety of contexts 
(see e.g! summaries of applications). It was shown 

in Refs. [45j HH and [61] that observables in a vast range 
of problems of ID non-equilibrium interacting fermions 
(and bosons) can be expressed in terms of Toeplitz de- 
terminants A A? = det {9i-j) x <ij< N - 



The behavior of determinants of such matrices becomes 
particularly non-trivial when the corresponding symbol 
(essentially the Fourier transform of gj-k, see Sec. IV B 
for more detail) has singular points known as Fisher- 
Hartwig singularities. In our case such singularities arise 
in view of discontinuities of the double-step distribution 
function /++. In Refs. [BTl and R)2l the large- A behavior 
of Toeplitz determinants with Fisher-Hartwig singulari- 
ties has been established analytically and verified numer- 
ically. These results ("generalized Fisher-Hartwig con- 
jecture") go beyond the "standard" Fisher-Hartwig con- 
jecture (proven in Ref. I63p as they contain not only the 
leading term but also subleading power-law contributions 
that have different oscillatory factors. We will see be- 
low that taking into account such contributions will be 
crucial for obtaining the oscillatory dependence of visi- 
bility of MZI on voltage. A further generalization was 
achieved recently in Ref . [52] where a broader class of sin- 
gular Fredholm determinants (determined by two symbol 
functions that show multiple singularities in energy and 
coordinate spaces, respectively) was explored and cor- 
responding asymptotics were found. Such determinants 
will arise below when we will invert the operator 2?**. 



2. Asymptotics of Toeplitz determinants 

For the benefit of the reader we summarize here the 
relevant results on Toeplitz determinants. 

A Toeplitz matrix gjk with j,k = 1, 2, . . . , A is defined 
by its symbol g(z) as follows: 



gjk = gj-k = 



.2tt 



(94) 



The determinant of such a matrix is called Toeplitz deter- 
minant. An important class of Toeplitz matrices (which 
is of relevance for our work and for various other non- 
equilibrium many-body problems) is generated by sym- 
bols with Fisher-Hartwig (FH) singularities, 

m 

g(z) = e v ^ ]J\z- Zj \ 2 ^ l3 (z) (z/zjf* (95) 
j=o 

where V(z) is a smooth function, m + 1 is a positive 
integer (number of singular points), Zj = e l<Pj , Reaj > 
— ^, pj € C, and 



— 7r < argz < ifj, 
ifj < arg z < tt. 



(96) 



In the context of our work, only the case a,j = (when 
the singularities of f(z) are discontinuities) will be rel- 
evant, so that we consider it henceforth. The large- A 
asymptotic behavio r of the corresponding Toeplitz de- 
terminant An read d 61 1 6 ^ : 
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no + ...+n m =0 j=0 



0<j<k<m j=0 



(1 + ...). 



(97) 



r 



/^^^(e^) and G is the Barnes G- 



where Vq 

function. The summation in Eq. ( |97[ ) goes over a set of 
integers no, n\, . . . , n m (whose sum is zero); we will see 
below that they can be understood as labeling branches 
of log g(z) in the intervals of continuity of the symbol. 
Each of these sets ( "branches" ) is characterized by a dis- 
tinct factor rij=o Z J 3 that m our con text will give rise 
to a distinct oscillatory exponent. Equation (97) presents 



explicitly the leading asymptotic behavior for each of the 
branches. There exist also subleading power-law correc- 
tions within each of the branches (i.e., corresponding to 
the same oscillatory exponent); they are abbreviated by 
+ . . . in the last bracket. Such corrections will be of no 
importance for our consideration, and we discard them 
below. 



We return now to determinants of the type ( 59 ) that 



arise in the course of the study of MZI. Here / is some 
distribution function and S(t) = 5wm !T ](t) is a constant in 
the window of the duration r and zero otherwise. We are 
interested in the large-r asymptotic behavior of A r [5, f] . 
As discussed in Sec. 



IVB1 



the UV regularization is im- 
plemented by using a high-energy cut-off A, so that the 
energy is restricted to the range [—A, A] and the time is 
discretized, tj = jAt = jn/A. In energy representation, 
the operator of interest reads [cf. Eq. (93)] 



To be specific, let us consider explicitly two examples 
corresponding to two lowest values m — 0, 1 (i.e., one and 
two FH singularities.) First, we consider the equilibrium 
distribution function /(e) = 8(fi — e). The symbol is 

a(e^) I ^ /(2 ^ 5 , -tt < <p < Tr/x/A, 



which is of the form (95) with m = 0, a?o — 0, /3q = 
5/(271-), z = e"^ /A , and V = iS(l + fi/A)/2. According 
to Eq. (97) in the large- N limit the det^-fc) asymptot- 
ically behaves as 



A [5, /single] = exp 



Z7T 



Ar 

7T 



xG(l-f )G(l + f ). 

Z7T ITT 



(101) 



Next, we consider a double-step distribution function 
/(e) = (1 — a) 9(po — e) + a6(ni — e), where we assumed 
that /io < fJ>x. In this case the symbol reads 



JS V /(2ir) 



A 
A 



A ' 
< (f < IT. 

(102) 



3(e) = 1 + (e" 5 - l)/(e). 



(98) 



Hence, the symbol has two FH singularities Zj — e 47r w/ A , 
j = 0, 1, with 



This can be identified with a symbol g{z) of a Toeplitz 
matrix, provided energy e G [—A, A] and angle ip € 
[— 7r, tt] are related by rescaling: ip = en /A. The intro- 
duction of a hard cutoff ±A and the above compactifica- 
tion of the energy axis will give rise to unphysical effects 
at this energy scale (since it will generate an unphysical 
discontinuity at ip — ±tt). These unphysical effect are 
eliminated by imposing "periodic boundary conditions" 
in energy domairPD which amounts to the following mod- 
ification of the symbol: lim e _ > ._Ag(e) = lim^A ff(e): 



7(e) =e 



We/(2A) 



[1 



JS 



l)/(<0] • 



(99) 



Here we have taken into account that lim e _i._A /(e) = 1 
and lim e _j.A /(e) = 0. 



-2irif} 



We choose 



^8 



l)a 



-2m0i 



1 + (e lS - l)a 

(103) 



/3i 



In 1 1 + (e iS - l)o| , fa 



2tt 



2tt 



Pi- (104) 



It is easy to see that the symbol has the form ( 95 ) with 
m = 1, otj = 0, and 

V(z)=V =iS/2 + iS^+ieVjPi (105) 

where we introduced eV = /ii — /io- According to 
( [97] ) the asymptotic behavior of the Toeplitz determinant 
det (gj-k) is given by 
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A [8, /double] = exp 



x pVt 
i_(A + / xo)r+— ln[l + (e«-l)a] 



E 



-ieVrn 



At 

7T 



-(/3o+n) 2 -(/3i-«) 2 



neV 



2(/3 +n)(/3i-n) 



G(l + /3 + n)G(l - A) - »)G(1 + ft - n)G(l - A + »). (106) 



In order to identify in the sum over n the leading contri- 
butions in the long-r regime, we consider the exponent 

Re [-(A) + n) 2 - (ft - n) 2 ] = -2 (n - n*) 2 + const, 

(107) 

where 



= -Re(ft - A>) 



= — Im In |7l — a) 
2w LV ' 



S_ 
47r' 



(108) 



Note also that the sum of voltage and time exponents, 
[-(A, + nf - (ft - nf] - [2(A) + n)(ft - n)] = ~(ft + 
ft) 2 is independent of n. Thus, terms dominant for 
At ^ 1 are also leading for large voltages, eVr » 1. 
For the analysis of the optimal value n* , we make the 
decomposition 8 = 2irM + 8' with MeZ and \8'\ < tt. 
One can show that the phase 



8" = Im In [(1 - a) + ae lS ] 



(109) 



has the same sign as 8' and satisfies \8"\ < \8'\. Then the 
optimal n* becomes 



M 8' - 28" 



11 = Y 



4tt 



1 7i* +M/2\ < 1/4. (110) 



We see that in the case of even M there is a single contri- 
bution with ri = —M/2 giving the most significant contri- 
bution to the asymptotic series; other contributions have 
substantially smaller (by real part) exponents. On the 
other hand, for odd M one has to take into account two 
contributions with n = —(A/ ± l)/2. Indeed, if a = 1/2 
(and thus 8" = 8' /2), then these two contributions come 
with equal exponents. When a deviates from 1/2, the 
exponents become different but still may be very close. 

It was shown in Ref. [52] that these results can be gen- 
ralized to a broader class of singular Fredholm determi- 
nants. Specifically, consider a matrix 



9j.k 



dc 
2A 



e -ieTz/A\y-k-S(t j )/(2w)]~^. 



with the symbol 



i 



l) /(e)- 



(Ill) 



(112) 



Here the notion of symbol has been generalized to in- 
clude both time and energy dependence (through the 



function 8(t) and /(e), respectively). Let us focus on 
the case when both the phase 8(t) and the distribu- 
tion function /(e) are piecewise constant functions with 
jumps at times t\ < t-i < . . . < tn t and energies 
Ml < A*2 < ••• < Mw,,, respectively. They satisfy the 
boundary conditions 8(t) = for t ^ [ti,tjv t ], /(e) = 1 
for e < /ii, and /(e) = for e > fijf . The UV cutoff and 
the periodic boundary conditions in energy domain can 
be implemented as before. The discontinuity points de- 
fine a grid which subdivides the time-energy plane in do- 
mains with different values of the symbol. The domains 
can be labeled by the time indices j £ {0, . . . , N T }, and 
energy indices k £ {0, ... , N^}. One associates with this 
set of domains a set of number Cjk, 



1 



Cjk 



lng(Tj + 0, n k + 0) + rijk, 



(113) 



2m 

Cjo = 8{tj + 0)/(27r), c ok = c Nrtk = Cj tNft = 0. (114) 

where {fijk} is an arbitrary set of integers. Further, a 
matrix with a time index j £ {1, . . . , A^ T } and energy 
index k £ {1, . . . , AL} is introduced according to 



Pjk = c j.k-l — Cj.k + Cj-i.k 



-i,fc-ii 



(115) 



Physically, each entry of this matrix corresponds to a 
crossing point of one energy-space and one time-space 
singularity. In terms of this matrix, a set of time (pji) 
and energy (qkm) exponents is defined as follows: 



Pjl = E Pjm'filr, 
rn' — 1 



N t 

E 



(116) 



One should keep in mind that Cjk and thus /3,'fci pg, and 



Qkm depend on the set of integers rijk- In Eq. (113) the 



logarithm ln^ is understood as evaluated at its principal 
branch, Im In g G (— 7r,7r]. The summation over integers 
rijk hence amounts to summing over different branches 
of the logarithms. 

We are now ready to state the result. For large time- 
and energy differences, \(tj - n)(^ k - Mm) I > 1 ( 3 h 
k ^=m), the asymptotic behavior of det(gj- k) is given 
b^ 
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A[S(t),/(e)]= r {n Jfe} exp 

{rijfc} 



l<3<N t 



l<k<N^ 

>< n n 

l<j<KN t l<k<m<N fi 



- Mfe) j {rj- 


HI - 


■Ti) 




A(r,-r ; ) 


Pjl 


tt(ij 


k f-^m ) 


7T 






A 



(117) 



with coefBcients T{ n . k y that are independent on Tj and 
Hk. It is not difficult to check that for the phase 8(t) 
proportional to a window function this formula agrees 



with the asymptotics (1061 of the Toeplitz determinant. 
While a rigorous mathematical proof of the asymptotic 
formula (117) is still lacking, Ref. 1521 presented powerful 



analytical arguments in favor of its validity supported by 
strong evidence based on numerical evaluation of such de- 
terminants. We will use Eq. ( |117[ ) below to get analytical 
results for the current through the MZI. 



3. Inversion of the single- channel counting matrix D„„ 



As we have shown in Sec. IV B the interference current 



can be expressed in terms of the inverse of the single- 
channel counting operator 2?**. We show here that it 
is related to a generalized Toeplitz determinant, whose 
asymptotic behavior can be estimated on the basis of 
results presented above. To this end we consider the 
time discretized expression for the operator I?**, given 



by Eq. (80), which has the symbol 

= 2iiJ + e/A 



9(e) 



1 



+ -1 /++(e) 



(118) 



corresponding to the Toeplitz matrix, 
i 1 



9j-k 



- (««. - 1) 



2tt j - k - 2t?+/ 
x " jRie -«*£>/Afa--k-2*+ M + Tl ~ 

The inverse of matrix g reads 

,det gi(k,j) 



det(g) 



(119) 



(120) 



where g^(k,j) is the (N — 1) x (JV — 1) matrix derived 
from g by removing the fc-th row and j-th column. 

Since our primary interest is t he matrix element 
D~^(t,t) with i < t < i+T, see Eq. (90), we concentrate 
specifically on the element (g~ 1 )ik, i-e. we put j = 1. In 
this case one has 



(s*(M)) 



lm 



gi. m +i, 1 < I < k, 

9l,m, k < I < N - 1 



(121) 

_ ( p m + {tv,t h ) _ A 
2„, [ e L ) 



k<l<N 
1 



x R rL e- inU / A ^>- k -^ tl ' t ^+Ti], (122) 



r 



where we have introduced the time-dependent phase 
,9 ft -t\ - I ^++^A t<U <t k , 



with t\ = t + (I — l)At. In the continuous representa- 
tion the phase $+(t, ife) is the piecewise function of time 
t. Tak ing i nto account Eqs. (119), ( 122| and the defi- 



nition (111) one observes that the matrix (g" (fc, 1)) lm is 



the generalized Toeplitz matrix with the symbol (112) 
where the phase 8(t) = 4$ + (i, i^). Its determinant can 
be dealt with by the use of results presented in the end 
of the previous subsection. Hence, 



(* _1 )i* = (-l) 



1+fc A[4fl + («;t fc ),/ ++ ] 
A[4tf + ,/++] 



(124) 



It is instructive to apply first the asymptotic rela- 
tions (106) and (117) to invert the matrix gj-k in the 
limit Ti -> 0, when the alternative evaluation can be 
done via the Riemann-Hilbcrt method. Under this con- 
dition the operator 2?** has the form 



J7W0 



= A, 



/o> + 



4id'' 



./if 



a: 



with A+(i) = 



Therefore [D a 



*JTi->0 



(125) 



can be evalu- 



ated in the similar fashion as we have found the inverse of 



the counting operator 25 in Sec. IV A 2 By introducing 
the function 



Y(t) = 



t-t 
t-i-r 



2i? + /tt 



(126) 
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-4n/v 



t 4 



In 



t+T 



FIG. 9. Phase 8(t) — 4i9 + (t;i fe ): on top of the "background 
phase" 4i?+(t) the 2-7r-pulse accounts for the additional elec- 
tron during the time interval < t < tk- 
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which solves the RH problem Y_ 1 Y + 
ther represents 2?** in the equivalent form 



-4ii) 



+ , one fur- 



*JTi->0 



= a + f_ 



Yl'ft+Yr 1 /?] AT 1 . (127) 



Using now the relation (85|), we obtain 



-i /+ _ r- 1 ^). 
/ j t,t 

(128) 

Taking into account the explicit form of the solution of 
the RH problem ( 126 1, for times t £ [i, t+r] we eventually 
arrive at 



r 1 



A 

•— e 

7T 



(t,i) = e - 2 ^+ sin(2<?+) e- iey( '~ 4) 



x It-tl-^+A-^f+r-tl^+A 

x T -2*+/* Ar »3+/*-i (129) 

where the regularization Yl (t) = Y_ (t+A^ 1 ) was chosen. 

Let us now demonstrate that the same result can be 
derived using the properties of the Toeplitz determinants. 
For T>~^(t,tk) and thus (g~) lk one needs to consider 
the generalized Toeplitz problem with 3 jumps in time 
domain, t\ = t, r 2 = tfc, T3 = i + r, and just 1 jump in 
energy domain [i\ = eV. Further it is cio = 2 r d + /n + 1 
and C20 = 2$ + / 7r an d hence pi 2 = —2d + /ir — 1, P13 = 



- (2tf+/ 7r + l) 2$ + /tt andp 2 3 = 2tf + /n. The asymp totic 
behavior of det(g) = A[4z9^,/ + ] is given by Eq. (101) 
with /i = eV and S — 2d + /ir. As the result, one obtains 



= Te 



-ieV(t-t k ) 



A(t k - t) 



-2t? + /7T-l 



A(t + T~t k ) 


2i5 + /tt 


At 


-2i9 + /tt 


7T 




7T 





(130) 

Except for the dimensionless unknown factor T and the 
prefactor A/ir = (At)" 1 which arises due to time dis- 
cretization, the above asymptotics agrees in all power- 
laws with the exact result (129). 



We now turn to the general situation of arbitrary T\. 
The distribution function in this case is given by / ++ 
instead of /+ which adds a discontinuity at \i\ = (the 
one at eV is now denoted by 112)- The asymptotics of 
the determinant A[4$ + (»; t k )] is determined by Eq. ( |117[ ) 
with 

Cio = ai + l, c 20 = «i, c kl =/3i+n k (fe = l,2), 

(131) 

where we have abbreviated 



ai = 2d+/Ti, ft 
and the exponents 



1 

27TI 



In 



(132) 
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-Pi- ni)(ni 


- n 2 


-1)4 


- (A + ni)(n 2 - ni), 


P23 


= "i - 


- n 2 


-l)(/?i+n 2 - 


Oil) 


- (n 2 


- nx)(px + n 2 ), 


Pl3 


= (H 


- OL\ 


-»i)(y8i 


f «2 


- «i) 


- (/3i +ni)(/3i +r»a), 


qi2 


= (H 


- Ct\ 


-/3 1 -n 1 )(/3 4 


- ^l) 


+ (ni 


- n 2 - l)(n 2 - ni) - (/?! - ai + n 2 )(/3i + n 2 ) 



One then has 

(-l) 1+fe A[4tf + (.;i fe ),/ ++ ] = r (ni ,„ 2) e MlAr+lAeyr e leV ' [ni(tfc -* )+Tl2(f+r - tfc)1 

ni,n 2 



A(t fc - 1) 


P12 


A(t + T-tfe) 


P23 


At 


P13 


ireV 


7T 




7T 




7T 




A 



(134) 



When deriving this asymptotics, we took into account 
that the phase factor e lA (tk-t) _ l) 1 + fe j since t k = 
i+ (k — 1) At with the infinitesimal time increment At = 



7r/A, cf. Eq. (123). The above relation is one of the main 



results of the section. It yields the a symptotic value for 
(g _1 )ife, expressed through Eq.(124). The determinant 



A[4$ + , /++] ap pearin g in the latter relation can be found 
exactly via Eq. (106), where one has to set 6 = Ad + and 
a = R\. 

The following remark must be made concerning the 



above calculations. The result ( 134[ ) has been derived 
with the use of the asymptotic formula (117). It is valid 
under the assumption \(Tj — Ti)([i k — /i m )| ^> 1, which 



defines the range of applicability to Eq. ( 134 ) , namely 
tk —t> 1/ eV and i+r — t k > 1/ eV . Below we examine 
another limit, when the time t k is close to either of two 
boundaries, t or t + t. 

To this end we represent the (normalized) determinant 
A[5(t)J(eM = A[6(t), f(e)]/A[S(t),T = 0] in the equiva- 
lent formal 
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A[<K*),/(e)]= J2 f {«,-} ex P 



{n 3 k} 



E E 

L i<i<JV t i<fe<Af^ 



jftfcMfc 



n n 



(135) 



where 



(136) 



The normalized determinant is cut-off (A) independent. 
All dependence on A comes from the zero temperature 
determinant, which up to a constant prefactor reads 



A[S,T = 0] = exp 



E At J 

l<j<N t 



2tt 



(137) 



>< n 

l<j<J<JV t 



A(Tj -7j) 



(5^-aj_i)(<5,-5,_i)/47T 2 



where we have defined the phases (5, = <$(£,■ +0) . Eq. ( 137 ) 
is the particular case of the asymptotics ( |117 ) when the 
generalized Toeplitz problem has N t jumps of the phase 
S(t) in the time domain and the single (Fermi) edge at 
e = 0. The summation over branches of logarithms is not 



required here. The equivalence between two forms of the 
asymptotic expansion, Eqs. (1371 and (117), follows from 
the sum rules, 



$>* = 



k=l 



2tt 



E^ = °- 

3=1 



(138) 



As before, representation (117) holds provided (t; — 
Tj)(l^m ~ Mfc) ^ 1 for all j > ^ and m > k. However, 
it enables a natural generalization to the situation, when 
this condition is not satisfied. Namely, if for some set 
(l,j,m,k) the opposite inequality is fulfilled, the corre- 
sponding factor has to be omitted from the product in 



Eq. (117). This is the advantage of the normalized rep- 



resentation in comparison to Eq. (135). In this way one 
can find the asymptotic form of (D~J-)(i, t) if t is close to 
t or t + t. 

If tk — i < 1/eV, we obtain 



(-l) 1+fe A[4tf + (.;t fe ),/ ++ ] =J2K 2 e laiAT+l(Pl+n2)eVT 



A(t fc - 1) 
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At 
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neV 
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7T 
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where we have introduced the exponents 

Pu = -(1 + ai)i Pis = -("i -Pi- n 2? - (ft + n 2 ) 2 , q' 12 = 2(ai - ft - n 2 )(ft + n 2 ). 
In the other limit, t + t — < 1/eV , the asymptotic expansion yields 

«" 2 



(-l) 1+fc A[4tf + (.; **)> /++] = E r:; i e MlA ^ i ^+" 1 



A(i + T- tfc) 


P23 


At 


Pl3 


ireV 


7T 




7T 




A 



with the exponents 

pgg = on, p'/g = -(1 + ai - ft - m) 2 - (ft + m) 2 , q'( 2 = 2(1 + oi - ft - m)(ft + n x ). 



(139) 



(140) 



(141) 



(142) 



We notice, that Eqs. ( 139 1, ( 141 1 and ( 134 ) represent the 



asymptotic expansion of the generalized Toeplitz deter- 
minant A[4$ + (»; tfc), /++] in the different domains of the 
variable tk- It is straightforward to check, that asymp- 
totic formulae (139) and (134) match each other at the 
scale tk — i ~ 1 JeV because of the mutual relations 

(143) 



Pl3 =P13 + P23, 



1X2 - Pl2 = 912 - P12 



between the power-law exponents. Similarly, the ex- 
pansion (141) matches Eq. (134) at the time scale 
(t + t — tk) ~ 1/eV due to analogous relations 

Pl3 = P13 + P12, q" 2 - P23 = 912 - P23- (144) 



The sketch of the determinant A[4$+(»; tk), /++] as the 
function of the time tk is shown in Fig. |10| 



To summarize this section, we have found the dis- 
cretized representation [g~ 1 )ik for the inverse of the 



single-channel counting operator X>**(t, tk), see Eq. (124). 
The generalized Toeplitz determinant A[4?9 + (»; tk), /++], 
appearing in this formula, is given by the asymptotic 
Eqs. (134) and ( 139|141 ). Accordingly, the denominator 
A [4$+, /++] can be found with the use of Eq. (106). 
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C. Interference current in the strong coupling limit form 



In Sec. |IVB| we have discussed the asymptotic prop- 
erties of singular Fredholm determinants and have found 
the inverse kernel V~^(t,t) of the single-channel counting 
operator. We are now going to use these results to eval- 
uate the AB conductance in the strong-coupling limit. 

We start by considering the particle number -ZV + _, 



which is given by Eq. ( 42 ) of Sec 
the relation -4f er m 



IIIB 



Making use of 
UnDeW together with Eq. Q, 



we can represent 



NH__l_ in the form 



N + ^=-i(R 2 T 2 )^ 2 , 



/<!> 



dt e 



DetX> 
DetX>(°) 



-\fV-\i,t)\-}, 
(145) 
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fen 



SA 



where we have defined the action A^ 

Ti=Q 

in the absence of edge- 



and the operator 



phases 



to-edge tunneling. Eq. ( 145 1 is evaluated at the optimal 
q (i) = d v Wt }T {t) with $ v = —rjir/is, which have 

It 



different sign for upper and lower arms of the MZI 
turns out that the integrand is in fact independent of 
time t and thus the integral is formally divergent. This 
amounts to an infinite number of electrons counted dur- 
ing an infinite measuring time in a stationary situation. 
The stationary current its obtained by dropping the i- 
integral and putting, say, t — 0. Using Eqs. (88), (90) 
one arrives at 



t + T 



DetX>(°) 

dt' 2?" 1 (i, *')/+- (t'-t) (146) 



Let us now use the result of Sec. IIV Al where the evalu- 
ation of the matrix determinant Det T> has been reduced 
to the product of single-channel determinants. For the 



scalar operator 2?** we introduce T> 



(0) 



V, 



Tak- 



ing into account the factorization formula (78) and the 



21=0 



block structure of the matrix 2?* , given by Eq. ( 79 1 , one 
can write 



Det2? _ Det2?„ _ A[4tf + ,/ ++ ] 
Detm ~ Det p(o) " A[4# + ,/+] 



(147) 



(to derive this relation we have made use of the fact that 
the operator T)q is Ti-independent). The dimensionful 
operator kernel T>~^(t\, t 2 ) and its discrctizcd dimension- 
less counterpart g" 1 are related by the energy factor 



W 



where the label "(0)" denotes the T x -+ limit. With 
the help of Eq. ( 124 ) for the matrix element g7j~ the 
expression for the current 1^ is then reduced to the 



e~ x I+- = (R 1 T 1 R 2 T 2 ) 1 ^ 2 e '*+2^+ e ^ (0) [ T dtk 

Jo 

„ (-l) 1+fe A[4tf + (.;t,),/ ++ ] 



x W 



A[4^+,/+] 



(f+(t k )-f-(t k )). 

(148) 



Here we have substituted the off-diagonal density ma- 
trix element /_| = i(R\Ti) 1 l 2 (f + — /_). In the formula 

above, one needs to perform further the integration over 
time t k and to evaluate the action A^ in the absence of 
tunneling between two edges of the MZI. These steps of 
calculations are discussed below. 



1. The time integral over tk 

To perform the time integration let us consider the in- 
dependent part in Eq. ( |148 ), 



J(t k ) = (-l) 1+k A[M + (.;t k )J ++ }(f + (t k ) - /_&)). 

(149) 

According to the asymptotic analysis of the previous 



Sec. IV B this integrand has the power-law singulari- 
ties which give the dominant contribution to the inte- 
gral (148) around t k ~ and t k ~ t, provided the real 



part of the corresponding power-law exponents is nega- 
tive. In general, for any time t k £ (0, r) the integrand is 
a superposition of powers-law terms, 



(150) 



where the exponent 7 = p\ 2 +P23 +ftn — <7i2 ensures the 
correct dimensionality (which is inverse time). As dis- 
cussed previously, the sum here runs over integers m, n 2 
or both depending on whether the time t k lies in the re- 
gion II, III or I, respectively (see Fig. 10 1 . For instance, in 





t+T 



(eVT 1 



(eVT 1 



FIG. 10. Sketch of the singular Fredholm determinant 
A[4#+(»; tk), /++] as fu nctio n of time tk- In region I, the 
asymptotic expansion ( 134 1 is valid. In regions II and 



III, power-law exponents are different, and A is given by 
Eqs. ( |139| l and ( |141| |, respectively. 
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the region I the function J , in accordance with Eq. ( 134 ), 
has the above scaling behavior with p\2 = pn ~ lj 
P23 = P23, P13 = P13 and <7i2 = q i2 - Equation ( |133| ) 
shows that by choosing \ni\ and | tt-2 | sufhciently large, 
Kepi2 and Rep23, respectively, can be easily made neg- 
ative. Therefore the integral over t k will be determined 
by a vicinity of the end points of the time interval (0, r) . 

First, let us examine the limit of short times t k *C r. 
One has to consider two asymptotic regions. For A -1 < 
tk ^ ( e ^0 (region II) we can use the short time expan- 
sion 



/+(t fc ) "/-(**) = 



2tt 



1 



el/ 
2tP 



(151) 



which gives the powers p X2 = p' 12) P23 = 0, pi 3 = p' 13 , 
and qi2 = q[ 2 + 1. Evaluating the ifc-integral over the 
region II for some given integer ri2, we find 



1/eV 



1/A 



dt k J(t k ) ~ (eyrfi3( e F/A)«- +1 -f'i3 



x (At, 



^12 + ! 



1/A 
1/eV 



- (eyT)^(el//A) 1+Q i +ai , (152) 
where we kept only the dominant contribution. Here 



oi\ = -2/v and, cf. Eq. (1401 



Pl3 



2 ! — ^ i ) 2 -f (153) 



The term which gives the leading contribution to the 
current from the region II is then found by maximizing 
Rep' 13 with respect to n 2 . The fact that the exponent 
A is independent on n 2 is not a coincidence. It encodes 
renormalization effects due to high-energy virtual exci- 
tations. In contrast, the arbitrary integers which encode 
different branches of In g are relevant for intermediate en- 
ergies < e < eV only, and thus do not affect the high 
energy scale A. 

Let us further look onto longer times, (eV)" 1 < t k <C 
r from region I, where we can approximate f+(t k ) — 
f-{tk) ~ 1/tfe- Hence the powers are p\ 2 — Pn — 1, 
P23 = P23, P13 = P13, and q 12 = q\2- By choosing some 
intermediate time scale (eV)^ 1 < t -C r as the upper 
cut-off, the integral reads 



1/eV 



dt k J(t k )~{eVT)^ + ^ 



x (eV/K) qi2 - pw - p23 {kt k y^ 



1/eV 



Under the assumption Repi2 < 0, the upper boundary t 



is irrelevant. Using relations (140) and (143), we obtain 



exactly the same asymptotics as in Eq. (152), 



1/eV 



dt k J(t k )~(eVr) p ™(eV/A) 



(154) 



We turn now to the analysis of the integral ( |148[ ) 
around the second singularity t k ~ r. Close to this end 
point one has f+(t k ) — f-(t k ) ~ r^ 1 . Following the same 
line of reasoning as above, we consider two asymptotic 
regions (I and III). The time integral over the region III 
for a given integer n\ yields (we recall that we consider 
the case v > 2) 



-(eV)-i 



dt k J(t k ) ~ (eVT) p "°- 1 (eV/A) 2+a * +2ai 



l+ai 



X (A(T-tfc)) 

T-(eV)- 

~ (eVr) p "»- 1 (eV/A) 1+0 ^ +ai , 



T-A" 



(155) 



The exponent p'/g — 1 can be read from the defini- 



tion (142). It is convenient to rewrite it in the form 



analogous to Eq. (153) 



p" 3 - 1 = -2 m 



ai + 1 - 20i 



(ai + 1) 2 



(156) 



which explicitly shows that Kep'{ 3 — 1 can be maximized 
with respect to n\. 

It remains to estimate the time integral when t k lies 
in the region I. Introducing as above some intermediate 
time scale t satisfying (eU)^ 1 « ! < t - (eU) _1 (which 
in case of Rep23 < — 1 will be irrelevant as an integral 
boundary) one obtains 



-(evy 



dt k J{t k )~{eVTy^- l {eV/PC) 



(157) 



In the following we are interested in the integers n 2 and 
n\ which maximize the exponents, Rep' 13 and Rep'/g — 1. 
To this end we write a,\ — M + m with M € Z and 
\ m \ < 1/2- Then for even M leading contributions come 
from n 2 = M/2 and m = (M + l)/2 ± 1/2, and for odd 
M they come from n 2 = M/2± 1 /2 and m = (M + 1)/2. 
Straightforward analysis shows that for integer filling 
fractions v > 2 in all optimal contributions we have 

Rep' 13 > Rep'/a " !• 

These observations lead us to the conclusion that, with 
all oscillatory terms e leVTn 2 e ia i Ar an( j ^-independent 
contributions, 



At 
n 



(158) 



taken into account, the leading terms of the t^-integral 
for v > 2 are 
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f T 1M dt k J(t k )/A[M+J + ] = A e -i«ieVT+i0ieVT( eV/A y +ai 
ll/A K 



(159) 



Ed 



ieVrn 2 



where T' n2 and are some unknown dimensionless con- 
stants. This expansion contains all terms in leading order 
of (eV/A) (also those subleading in (eVr)). 



2. Action A^ in the absence of tunneling 

Let us now evaluate the action of the system when 
inter-edge tunneling is absent, 



iA (a) = TrLn 1 



f + e 2M "f -2iTn? 9 / . 



(160) 



In this subsection the traces extend over all v upper and 
v lower inner channels. We combined all 2v distribution 
functions f\ and phases d\ into 2v x 2^-matrices / and 
•d q . Due to the Dzyaloshinskii-Larkin theorem we antici- 
pate that only first-and-second-order-in-$ terms are non- 
vanishing for the action above (it is worth reminding here 
that throughout this section the distribution functions 
were assumed to be the Fermi-like). Hence we expand, 



Let us assume that among the channels belonging to the 
upper edge of the MZI, the outer channel is biased by 
V and all the rest ones by Vq. At the same time, all 
channels on the lower edge are grounded. This gives the 
first, "zero mode" contribution to the action, 



2iTr§ q f - 



(;-/o) 



.e[V + (v - 1)V ] 
1 2* 



2§ 4 



(164) 



The quadratic contribution to the action iA^ is UV di- 
vergent and needs to be cutoff by at the scale A ~ oj c , 
yielding 



-2Tn?« 1 - / d q f = 



2vd\ 



In uj c t 



lnc 



(165) 



In passing we note, that since is a purely Gaussian 
contribution, it can be equivalently evaluated by averag- 
ing the phase e tAm = (V<(o)-« b -(o)\ over the Gaus _ 

sian action of the MZI in the limit T x — > 0. 



iA {a) = Tr ^Ln ^1 + [2iti q - 2^ 2 j /J - 2i® q f 

= 2iTrtf« (J - / ) - 2Tr d q (t - f)S q f. (161) 

Consider now three local operators A, B, C, where by 
definition A(ti,t2) = A(t\)5(ti — t 2 ) etc. Evaluating 
the following trace (one should carefully take into ac- 
count here the non-local in time structure of the Fermi- 
distribution function), one obtains 



Tr t A[BJ ]C = fdt 2 lim A(tj) 



i B(h)~B{t 2 ) 
2tt h-t 2 + i0 



2n 



dtA(t)B(t)C(t). 



C(t 2 



(162) 



For Fermi-like non-equilibrium distribution functions, 



hit) 



-ieV t t 



fo{t), the above relation implies 



TvtAif, - f )C 



eVi 
2tT 



dtA(t)C(t). 



(163) 



3. Current, Conductivity and Visibility 

The results of previous subsections enable us to evalu- 
ate the Aharonov-Bohm current in the MZI. Setting the 
voltage of all outer channels to zero, V = 0, and defining 
the exponents 



p'{n 2 ) = -2 [n 2 



ai - 2/3! 



1 



(IT 



2 ) 2 

2 



at, (166) 



/ g 1 + l-2^ V 1 af 
p(m) = -2lni I ~2 + Y' ( 16? ) 

one can represent the coherent current contribution in 
the form 

J coh - 2Rc/+_ = ^(i? 1 T 1 i? 2 T 2 ) 1 / 2 2Ree 4 * T , (168) 

where Iq is the normalized amplitude of the Aharonov- 
Bohm current. Collecting Eqs. (|146[), (|159|, (|164|) and 



( 165 1 together, we finally arrive at 



In 



2e 



-2iri/v e i(p 1 + l/u)eVT 



EC 



(eVr) 



p'{n 2 



Ed' 



-ieVr 



1) (eVr) 



p"(ni 



(169) 
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The following table shows the power-law exponents cor- 
responding to the terms which give the dominant contri- 
bution to the series ( 169 1 at each filling factor v: 



2 
3 

> 4 



leading powers 



p'(0)=p'(-l)=0=p"(0) 
p'(-l) and p'(0) 

p'(0) and p"(0) if Ti < 1/2 or p"(l) if T 1 >l/2 

Taking these leading terms into account we obtain the 
results presented in Sect. II B namely Eq. Q with the 
exponents ( 12 1 and ( 13 ). We have also checked the valid- 



ity of analytical asymptotics ( 169 1 by straightforward nu- 
merical evaluation of Eq. (1481 for the AB current. The 



perfect agreement between the analytical and numerical 
approaches, demonstrated in Fig. |4j provides additional 
support towards the conjecture of Ref. [52] 

We close this section by providing the reader with a 
qualitative physical picture, underlying the result (1691. 
We begin with inspection of the phase pulse 5(F) = 
A r d + (t;tk) which according to our discussion in Sec. IV B 
determines the resolvent (2? -1 ) of the counting oper- 
ator (see Fig. 9|. First, we note that in our model 
of the maximally long-ranged Coulomb interaction the 
counting phase 4#+(t) = (Air / v)q(t) , with q(t) being the 
"background" charge on the upper arm of the MZI. The 
phase S(t) is different from 4i?+(t) by the 27r-pulse of the 
duration tf., which describes injection of the interfering 
"transport" electron at time i = into the interferom- 
eter via the QPC1 and its annihilation at time tk (the 
electrons leaves the MZI through the QPC2). As we see 
it from the calculations of the Sec. |IV C 1| in the high- 
voltage limit eVr 1, the critical exponents p' are asso- 
ciated with many-particle scattering processes with short 
tk, tk — h/eV. We can interpret it as the event when one 
electron enters the MZI through the QPC1 and another 
electron leaves the MZI shortly afterward on a time scale 
~ h/eV via the QPC2. The possibility for an electron to 
propagate through the system in a time tf. <C r is due 
to the long-range nature of Coulomb interaction in our 
model (a similar situation will occur in a model with a 
strong short-ranged interaction, where such a process can 
be mediated by the exchange of plasmons, which have a 
velocity exceeding by far the bare velocity of electrons). 
On the contrary, the second critical exponent p" is due 
to many-body scattering events with time tk — t. In this 
case the MZI is excited into a state with an extra elec- 
tron for a long time r » K/eV (it will be, e.g., the only 
possibility in the limit of weak short-ranged interaction 
when the behavior of the MZI is very close to the one 
with non-interacting electrons). 

Next, we associate the integers n\ and ni with the 
number of inelastically excited electron-hole pairs on the 
MZI in the corresponding time intervals (at < t < tk 
and t k < t < r, respectively). More precisely, one can 
interpret each term in the asymptotic expansion for the 
current J(tk) as the product of forward (electron-like) 
and backward (hole-like) many-particle scattering ampli- 
tudes, A b ({n~};tk) x A? ({n + };tk). An electron, propa- 



gating through the upper arm of the MZI leaves a trace 
in the bath of particle-hole excitations, which is encoded 
in the numbers and n\ of the excited electron-hole 
pairs after the electron's injection and annihilation. The 
dominant e-h pairs correspond to excitations from one 
Fermi edge to the other, so that the typical energy of 
an electron-hole pair is Yujj ~ eV. The corresponding 
many-particle amplitude should behave as 



A'({n+};i*) 



-inTeVtk— ITI2 eV(r — t&) 



Similarly, the backward amplitude is characterized by 
numbers and ■ In the case when the relative num- 



bers n\ = n{ — n 1 and n-i 



-n between the forward 



and backward scattering amplitudes are non-zero one ob- 
tains an interference term in the current with a phase 
which is linear in voltage. The two dominant terms with 
inequivalent phases in the series for the current J_j pro- 
duce the lobe pattern in the visibility of the Aharonov- 
Bohm oscillations. For example, at v = 3 such two terms 
are those corresponding to scattering processes charac- 
terized by the short time tk ~ h/eV but having different 
numbers ri2 = —1 and ni — 0, resp. 



V. NUMERICAL APPROACH 

In this section we present a numerical evaluation 
scheme for the interference current. Results, obtained 
within this scheme, corroborate and complement the an- 
alytical study of Sec. |IV| based on the particular form of 
the Fredholm counting operator T>. After making a few 
introductory remarks regarding the numerical formula for 
the Aharonov-Bohm current, which is suitable for a prac- 
tical implementation, we consider two cases when we did 
not succeed to obtain the analytical asymptotics and had 
to evaluate the Fredholm determinants numerically. This 
is, first, the case of non-equilibrium incoming distribu- 
tion (Sec. VA) and, second, the regime of intermediate 
interaction strength E c t ~ 1 (Sec. VB). 



As argued in Sec. |IV C[ the interference current is ob- 
tained from Eq. ( |42[ ) by dropping the divergent ^-integral 
and putting i = 0. For the purpose of the present section 
we rewrite this formula in the equivalent form 



J+_ = I + 



(eq) 



DetV 



fv- 1 (0,0)|-) 



(170) 



Here, 



(eq) 



-i {R2T2 



n!/2 i$ iAf eIIa +i5A 



\ff>-\o,o)\ 



eV=t Q 



denotes the near-to-equilibrium current, evaluated at 
very small voltages 6qt -C 1, while the label "norm" 
refers to quantities which are normalized with respect 
to near-to-equilibrium values, i.e. 



DetV 



= DetZ>/DetX> 



2G 



The density matrix / can be removed from the matrix 



element appearing in Eq. (1701 by using a relation 



i-xr 



e ***_l )ff)-\ 



which leads to 



fv-\o,o) 



<+| O-i(0, 0)R 



\v-i(p,o) 



eV=e 



By solving an appropriate Riemann-Hilbert problem we 
were able to show that close to equilibrium all effects 
of non-equilibrium and interaction, like dephasing and 
renormalization, are absent for the Aharonov-Bohm cur- 
rent. In agreement with the Landauer-Biittiker result it 
depends linearly on voltage, thus 



(°q) 



{R 1 T 1 R 2 T 2 ) 



1/2 a® eeo 
2tt' 



Corrections to this equilibrium expression are encoded in 
the normalized terms in Eq. (170 1. By discretizing the 



integral operator T> the latter become amenable to nu- 
merical evaluation. Below we discuss two different dis- 
cretization schemes. In practice, the choice between two 
depends mainly on the strength of Coulomb interaction. 



A. Non-equilibrium incoming distribution 

We consider the situation when the incoming elec- 
tron beam is diluted and is driven out of equilibrium 



even before scattering at the first QPC. More specifi- 
cally, we focus on the setup (that has been realized ex- 
perimentally) with an additional QPCO with a trans- 
parency < R < 1 placed outside of the interferom- 
eter and diluting the incoming beam. The incoming 
then distribution function acquires the double-step form, 
/+(e)= T 9(-e) + R a 9(eV-e). 

As in Sec. |IV| we focus here on the limit of strong 
interaction E c t ^> 1, where the counting phase (54 1 is 



a window function and the counting operator T> is of 
block Toeplitz form with a two-channel structure. Note 
that its decomposition into single-channel operators that 
was performed in Sec. |IV A"| requires trivial incoming dis- 
tribution functions / + , /_, i.e. zero-temperature Fermi 
distributions, possibly with different chemical potentials. 
Thus the non-equilibrium form of f+ constitutes an ob- 
stacle for a further analytical evaluation. At this stage, 
we do not know whether there is an analytical way to 
overcome this problem. We thus resort to a numerical 
evaluation of the determinants. 

To proceed numerically we use the same discretization 
procedure as in Sec. |IV B| where the energy cutoff A is 
introduced and the symbol g(e) of T> is required to satisfy 
periodic boundary conditions in energy domain, 



t + (e 2M - l)/(e) . 

Now the symbol g{t) £ C 2x2 has an additional 2x2- 
channel structure, the same holds for the diagonal ma- 
trix •& with diagonal entries In the discretized 
representation the counting operator T> then reads 



A 



de 



-itf [j-k] 



A 2A 

e 2zS - 1 



9(e) 



2tt 



j-k- j?/tt \-i {Rinf 2 R (e~ iv ' 



']+T + T 1 E i{R 1 T 1 ) l/2 R (e- i ^l 



j-fc+i9 + /7r 



TiR e 



[j-*+*+/*J +Ri+TiTq 



(171) 



The resulting visibility which follows from the above 
approach is shown and discussed in Sec. II B see Fig. [5] 
A good convergence was already achieved for moderate 
matrix sizes with N = At/it ~ 100. While we do not 
have a complete analytical form of the Aharonov-Bohm 
current in this case, the dephasing rate can be deduced 
from the leading large-r asymptotic behavior of Det T> ~ 
Making use of the results for the block Toeplitz 



-r/ r * 



determinant^', we arrive at 



de 
2^ 



Re In det 



i. /• -;<-'"-,/ (172) 



eV 
2n 



Re In 



T + Rq [R 1 e~ i2 ' n,v + T 1 e i27r /" 



In the case v = 2 the dephasing rate is simplified to = 
— (eV/2ir) In |2i? — 1|, as has been already mentioned in 
Sec. 1111 
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B. Intermediate interaction strength 



We now consider the discretization scheme applicable 
in the case of a moderate charging energy, E c t ~ 1. In 
this general case the correlation function (511 has to be 
evaluated numerically (see, e.g., Fig.|8]in Sec. IIIB which 
shows the result for lu c t — 25). Owing to finite E c , the 
counting phase $ + (t) is not a piecewise constant function 
of time anymore, but rather acquires oscillations in time. 
As an illustration, the phase 9+{t) is plotted in Fig. [7] 
for two different charging energies. Thus, the situation is 
different from the previous subsection, since the count- 
ing operator is no longer of Toeplitz form. Despite this 



is nevertheless possible, but a time discretization of the 
kernel T>(t,t') has to be performed directly in the time 
domain without any reference to the conjugate energy 
representation. To this end, we rely on the approach 
similar to the one used in Ref. (see Supplemental Ma- 
terial of that work) and interpret the zero-temperature 
Fermi distribution function / m terms of the Cauchy 
principal value and the Dirac delta distribution, 



fa{t) = 



2nt 



J- = ±*Mf(*). 
+ iO 2tt t 2 w 



complication, a numerical treatment based on Eq. (170) with tj = (j — l)Ai and At = 7r/A yields 



Consequently, the discretization of the non-equilibrium 
single-particle density matrix, fjk = sxf l (U,tj)s\At, 



fjk = (1 - Sjk) 



i 1 1 

1 04k 

2tt j - k 2 J 



Ri 
-i(RiTi) 



1/2 



CRiTi) 



1/2^ 



(1 - Sjk) 



i e 



2tt 



j-k 



Sjk 



eV 
2A 



r 



with A S> E c being a high-energy cut-off. Making use 
of above expression one can further construct the dis- 
cretized matrix of the counting operator 



V 



s jk t 2 



i) ° fjk 



(173) 



where "o" denotes matrix multiplication with respect to 
channel indices. 

The discretized form outlined above is very general, 
since it allows for arbitrary time-dependent phases. Note, 
however, that for the case of window-function time de- 
pendence of the "counting" phase (i.e., Toeplitz case) this 
regularization yields a result which is manifestly 2ir peri- 
odic in 2{t q , with nonanalyticity at points (2n + l)7r. This 
should be contrasted to the analytic and non-periodic 
behavior of the Toeplitz determinant within the proper 
regularization discussed above. The corresponding dif- 
ference between the present problem and that of FCS 
has already been mentioned in the end of Sec. IV A 1 
We have checked that for |2i9 9 | < ir both regularization 
schemes, Eqs. (171) and (173) produce identical results, 



although the convergence of the second scheme is gener- 
ally worse. 

In the case of a finite charging energy the phase 
changes continuously with time, and there is no more 
problem with the regularization (173), independently of 



how large the values acquired by the phase are. We used 
the matrix size N = At/it ~ 500, which was sufficient 
to obtain numerical results for the visibility at E c ~ 1 /r 
with good precision. The results were presented and dis- 
cussed in Sec. Ill CI 

Finally, we have calculated analytically the long-r 
asymptotics of the action Af crm = — ilnDetP, which 



gives the out-of-equilibrium dephasing rate (17) of the 
Aharonov-Bohm oscillations. 



VI. SUMMARY 

In this paper, we have discussed an exactly solvable 
model of a quantum Hall electronic Mach-Zehnder inter- 
ferometer for arbitrary integer filling factor v. The model 
is specified by a form of e-e interaction restricted to the 
inner part of the interferometer and two single-particle 
scattering matrices of quantum point contacts (QPCs). 
Our main results can be summarized as follows: 

(i) Making use of the non-equilibrium functional 
bosonization approach, we have established the ex- 
act solution of the above model in terms of the re- 
solvent of the Fredholm integral operator — single- 
particle "counting operator" T>, which is related 
to the problem of electron full counting statistics 
(FCS). The time-dependent scattering phase i?+(t) 
of the operator T> encodes all information about the 
interaction in the system. The link between the 
initial many-body problem with Coulomb interac- 
tion and single-particle quantities is established by 
virtue of the real-time instanton technique which 
becomes exact for the specific type of the Keldysh 
action describing the MZI. 

(ii) The focus of our study was on the model with "max- 
imally long-range" Coulomb interaction character- 
ized by the electrostatic charging energy E c . In the 
limit of strong interaction E c 3> 1/r (here r is the 
electron flight time through the MZI) the scattering 
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phase z9+ (t) becomes a piecewise constant "window" 
function and the operator T> simplifies to the block 
Tocplitz form. In the absence of external dephas- 
ing, we were able to get rid of the matrix structure 
of V and have expressed the result in terms of sin- 
gular Fredholm determinants that may be viewed 
as a generalization of Toeplitz determinants with 
Fisher-Hartwig singularities. This has allowed us 
to evaluate the Aharonov-Bohm (AB) conductance 
in a closed analytical form. At a moderate charging 
energy E c ~ 1/r and/or in the situation when the 
incoming distribution is made non-equilibrium by 
an additional QPC placed outside of MZI, we have 
obtained the results for the visibility by evaluating 
the determinants numerically. 

(iii) Results of our theory at E c <~ 1 /r match in all prin- 
cipal aspects the experimental observation in many 
designs of Mach-Zchndcr interferometers at filling 
factor v—2. If the transmission coefficient T\ of 
the QPC1 is close to 1/2 the visibility dependence 
on external bias shows a number of "lobes", their 
amplitude is being suppressed with the increase of 
voltage. The AB-phase dependence is close to a 
piecewise constant function with jumps equal to n 
at minima of the visibility. The visibility is further 
suppressed when the MZI is subjected to an out- 
of-cquilibrium shot noise, generated by the QPCO 
placed outside the interferometer. We have quanti- 
fied the dephasing rate which governs the de- 
cay of AB oscillations with bias in terms the trans- 
mission of QPCs, filling factor v and the strength 
of e-e interaction E c t. 

(iv) Our analytical results in the limit of strong interac- 
tion E c t >• 1 show an intimate connection between 
the observed "lobe" structure in the visibility, on 



one hand, and multiple branches in the asymptotics 
of singular integral determinants, on the other hand. 
In more physical terms, this is the many-body in- 
terference effect resulting from the quantum su- 
perposition of many-particle scattering amplitudes 
with the mutual phase differences which are linear 
in external bias. We derived the non-equilibrium 
quantum critical exponents, which depend both on 
the transmission T\ and the filling factor v. They 
are attributed to the Anderson orthogonality catas- 
trophe under out-of-equilibrium conditions and de- 
scribe the power-law dependence of the above many- 
particle amplitudes on voltage. 

Before closing, we mention some future research direc- 
tions. An extension of the presented approach to the frac- 
tional quantum Hall edge states devices, comprising two 
(or more) QPCs which couple the co-propagating edge 
modes, would be of great interest. Another important re- 
search direction is the analysis of the asymptotic behavior 
of block Toeplitz determinants with Fisher-Hartwig sin- 
gularities (and, more generally, block determinants, with 
symbols that have multiple energy and time singulari- 
ties). This would not only permit to obtain closed an- 
alytical results in the model with moderate interaction 
strength but would likely have multiple further applica- 
tions. 
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